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ABSTRACT 


The  non-linear  electromagnetic  theory  proposed  by 
Dr.  Schiff  is  investigated.  Nbn -Maxwellian  plane  wave  solutions 
of  the  field  equations,  k  =  -  (l  +  A^)  A^  ,  are  obtained. 

Static,  particlelike  solutions  of  the  field  equations  are  considered 
in  the  absence  of  a  magnetic  field.  Existence  of  well  behaved 
solutions,  in  the  case  of  spherical  symmetry,  is  formally 
established.  The  list  of  known  spherically  symmetric,  neutral 
and  charged  particlelike  solutions  is  extended  to  include 
"compound"  solutions  which  correspond  to  an  imaginary  vector 
potential  in  certain  (linear)  regions  of  space.  Sufficient 
conditions  are  derived  under  which  variational  approximations 
yield  upper  bounds  to  the  Lagrangian  associated  with  the  class 
of  field  equations  Ad  =  F’ (d) .  The  usual  parameter  variation 
method  is  generalised,  in  the  case  of  more  than  one  dimension, 
to  admit  variation  of  functions  of  one  of  the  independent 
variables.  This  method  is  used  to  obtain  approximations  to 
possible  nonspherical,  odd  parity,  neutral  particlelike 
eigensolutions  of  Ad  =  d  -  d^.  The  results,  however,  seem  to 
give  an  approximation  to  the  lowest  two-particle  state  rather 
than  a  one -particle  state  of  odd  parity.  A  variational  principle 
for  charged  particlelike  solutions  is  developed  and  used  to 
obtain  approximations  to  the  lowest  (in  energy)  spherically 
symmetric  state.  A  scheme  is  suggested  whereby  integral  relations, 
which  are  satisfied  by  integrable  solutions  of  Ad  =  FT(d),  may 
be  generated.  The  first  few  such  integral  relations  are  derived 
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and  are  used,  in  the  particular  case  F'(<t>)  -  4>  -  t J ,  to 
obtain  alternate  expressions  for  the  energy,  /  [ ^-(  V<t> ) " +  F(<t>)  ]d 
of  the  system  and  to  test  variational  approximations  to  the 
solutions  of  the  field  equation.  The  dynamical  stability  of 
neutral  particlelike  solutions,  against  certain  dissociation 
modes,  is  demonstrated;  in  particular,  against  direct 
dissociation  into  plane  waves  and  into  particles  for  which 
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I.  INTRODUCTION 


One  of  the  basic  problems  which  arises  before  every 
generation  of  physicists  is  that  of  "elementary  particles." 

The  physicist  attempts  to  determine  what  constitutes  the 
basic  "building  blocks"  of  matter,  how  they  arise  and  what 
are  their  properties. 

With  the  discovery  in  recent  years  of  a  great  number 
of  what  are  presently  considered  to  be  "elementary  particles" 
and  their  resonances,  this  problem  is  very  much  in  the  fore¬ 
front  of  research  both  experimentally  and  theoretically.  On 
the  theoretical  side,  quantum  field  theory  provides  a  good 
tool  for  describing  the  nature  of  the  interactions  between 
particles 

Recently,  the  investigation  of  symmetries  among  the 
various  properties  of  the  particles  (charge,  spin,  parity, 
isotopic  spin,  etc.)  has  led  to  an  apparently  successful 

3 

method  of  classifying  these  particles  in  symmetry  groups. 


1.  A  reasonably  up-to-date  list  of  observed  particles  and 
their  properties  is  given  by  A.M.  Baldin  and  A. A.  Komar 
in  Sov.  Phys . :  Uspekhi,  8,  p.  307^  Sep.  1965* 

2.  There  is  a  large  number  of  good  texts  on  this  subject, 
among  them  references  (l),  (2),  (3)  of  Bibliography. 

3.  References  (4),  (5)*  and  (6),  give  a  good  account  of  the 
present  state  of  this  rapidly  changing  subject. 
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This  in  turn  has  led  to  the  prediction  of  new,  hitherto 
unobserved  particles  such  as  the  ,  as  well  as  new  informa¬ 
tion  about  the  properties  of  "old"  ones.  Both  these  approaches, 
however,  do  not  lend  themselves  to  explaining  the  intrinsic 
nature  of  elementary  particles,  how  they  arise  and  why  with 
the  particular  values  of  such  basic  parameters  characterising 
them  (such  as  mass  and  charge)  as  are  observed  in  nature. 

One  approach  which  hopefully  can  shed  some  light  on  this 
problem  forms  the  subject  of  this  thesis.  We  consider  here 
the  classical  non-linear  electromagnetic  theory  proposed  by 
Dr.  Schiff  (7).  In  this  theory  the  basic  entity  is  the 
electromagnetic  field,  while  particles  (spinless  bosons)  are 
manifestations  of  this  field  arising  through  the  mechanism 
of  a  self  interaction  term  in  the  Lagrangian  and  hence  the 
field  equations. 

This  approach  dates  back  to  the  attempts  of  Mie^  to 
modify  Maxwell ls  electromagnetic  field  equations  in  a  way 
to  account  for  the  existence  of  particles.  The  hope  was  that 
the  modified  field  equations  w'ould  yield  solutions  which  differ 
from  those  of  Maxwell Ts  equations  principally  in  a  localised 
region  of  space  corresponding  to  the  "interior"  of  the  particle. 
One  of  the  most  widely  investigated  examples  of  Mie’s  theory, 
based  on  the  Lagrangian 


1. 


G-.  Mie,  Ann.  Physik  37.,  p.  511  (1912);  39,  p.  1  (1912); 

40  p.  1,  (1913).  A  good  account  of  this  theory  and  other 
attempts  is  given  by  Teshima  (8) . 
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was  considered  by  Born  and  Inf eld. 


1 


The  immediate  object  of  investigation  in  such  theor¬ 
ies  is  the  determination  of  particle-like  solutions  of  the 
field  equations  -  solutions  for  which  the  parameters  of 

the  field  and  all  physically  observable  quantities  derivable 

2 

from  them  are  well  behaved.  The  purpose  of  this  thesis  is 
to  examine  further  the  theory  proposed  by  Dr.  Schiff  (7). 

An  appreciable  amount  of  work  has  already  been  done  on  this 
subject^,  and  the  more  important  aspects  of  it  will  be 
abstracted  throughout  this  thesis  as  required  for  completeness. 

The  proposed  field  equations  are  (in  the  usual 
four-notation) : 


where 


and 


Fik 

An_ 


x 


k 


Ak,i  Ai,k 

(X,  it)  , 

(xj  ict )  . 


3 


1.  M.  Born,  Proc.  Roy.  Soc.  (London),  Al43,  p.  410.  (1934); 

M.  Born  and  L.  Infeld,  ibid.,  144. ,  p5  425  (1934),  147 j  P«  522, 
(1934);  150,  p.  41  (1935). 

2.  Throughout  this  thesis,  "good  behaviour"  will  imply  existence, 
continuity  and  single  valuedness  throughout  the  region 

under  consideration. 

3.  References  {j)3  (8),  (9 )f  (10) . 
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A  and  <D  are  the  usual  vector  and  scalar  potentials  respec¬ 
tively.,  while  m  and  g  are  universal  constants  having 
dimensions  of  mass  and  charge  respectively.  The  summation 
convention  on  repeated  indices  will  be  used  throughout,  with 
i,  jj  k,  ...  =  1,  2,  4  while  a,  (3,  7,  .  .  .  =  1,  2,  3* 


The  right  hand  side  of  (I-l)  is  identified  as  the 

electric  4-current,  which,  because  of  the  antisymmetry  of 

F. ,  ,  is  clearly  conserved, 
lk 

2 

me 

Under  the  transformation  -*■  •  ::3^—  and 

x  —x  •  —  ,  the  field  equations  take  on  the  convenient 
i  i  me 

(dimensionless)  form: 


F..  .  =  -  (l  +  A.  )  A, 
ik,k  v  k  '  1 


(i-la) 


In  vector  notation,  since 
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the  field  equations  are: 


-  J 


-o1.  al js JL  roq  •iati.-uH  baa  lotfoev  JE^t/au  erlz*  artjs  $  Los ’  -A 

■ 

<  .  ■ 


•  £-■*-1:  *1  “  -  ■  •  >3  siiriw  A  t£  tS  tI  «.  ...  <>i  «t  v-t 


-:  ...  f'x .  <  ^  .;• 


ot!  -  sol  s£  (l-l)  •  o,9i.>:  .  i  ;} •••.  *  *  • '  o.'P 


' 


Xl 


,A 


.  ■*.  ; .;  it  c  i  erij  ■  '/on  J 


y 


\  tc  9:  >&-••>'  ’  ■■  r  pi  frXdi'i  ortt  ,  - 

o/n 

^  jmiol  i  33 sX no:  andffllij ) 


-"(v 

+  X) 

• 

.  I 


•  ^ 


•:.t 


y r 


5 


Vxff  =  ||  -  (1  +  I2  -  02)  t 

(1-lb) 

V*E  =  -(1  +  T?  -  't’2)* 


The  theory  (as  is  obvious  from  (I -l))  is  relativistically 
covariant,  but  is  not  gauge  invariant.  The  conservation  of 
the  4-current  imposes  a  gauge  condition  on  the  4-potential 


(1  +  Ak2)  A±]  =  0  , 


which  will  be  referred  to  as  the  "generalised  Lorentz 
gauge."  In  vector  notation: 


(A •  V  +  0-^r)  ( 1  +  A 


-  *2) 


+  (1  +  A2  -  02)  (Vi  + 


^0  \ 

3t  ' )  " 


0 


Any  solution  of  the  field  equations  (i-la)  will  necessarily 
satisfy  this  gauge  condition.  We  note,  however,  that  under 
the  gauge  A,  2  =  -1  one  returns  to  the  Maxwell  free  field 
equations.  Except  for  this  case,  the  field  equations  are  non¬ 
linear,  so  the  principle  of  superposition  is  not  generally 
applicable.  As  a  corollary  this  implies  that  any  (stable) 
particle  will,  in  general,  decay  in  a  constant  potential  field. 
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The  field  equations  (i-l)  are  equivalent  to  the 
principle  of  least  action 


SS  =  0  , 


with  the  action  function  defined  as 


S 


i 

3 


where  the  Lagrangian  density  is 


-  h  ^  +  ?Af +  tW  + 


1 


(1-2) 


K  and  N  being  constants. 

The  symmetrised  energy-momentum  tensor  is 


-  If  c  +  W\2  +  K> 1 2 * *  5ik 

"  Fl£  Fk^  "  ^  +  AJl  ^AiAk^  *  (I_3) 


It  satisfies  the  usual  conservation  law 


T  =  o 

ik,k 


(1-4) 


1  Integration  over  the  entire  range  of  values  of  the  independent 
variables  x.  is  implied  throughout,  unless  otherwise  speci¬ 
fied.  d^x  =  dx1  dx2  dx^  dx^  throughout. 

2  Components  of  the  symmetrised  energy -momentum  tensor  are  listed 

in  appendix  A.  (1-3)  is  obtainable,  to  within. a_ divergenceless 

term  of  the  form  (A.K  «)«j  from  the  usual  definition 
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The  physical  implication  of  these  (four)  equations  is  that 
the  four-momentum 


j 


) 


where  2  is  any  spacelike  hypersurface  of  infinite  extension,, 
is  conserved.  =  ("?*  iW)  and  "?  is  the  total  (relativis¬ 

tic)  momentum  and  W  the  energy  of  the  system. 


Since  the  largest  transformation  group  under  which  the 
field  equations  ( I -1 )  are  invariant  is  the  ten  parameter 
inhomogeneous  Lorentz  group1,  then  according  to  Noether’s 
theorem^  there  should  be  a  total  of  ten  conservation  laws . 

The  four  conservation  laws  (1-4)  are  a  consequence  of  invariance 
under  four-translations.  Invariance  under  four  rotations  (the 
homogeneous  Lorentz  group)  leads  to  the  six  conservation  laws 


Mkim, k 


0  , 


(1-5) 


where 


M 


kirn 


=  (T,  -  T,  x„)  +  ( 

v  ki  m 


b£ 


baC 


km  $ 


31 


i,k 


m 


dA 


m,  k 


V 


+  s, 


(1-6) 


1  The  group  of  translations  and  rotations  in  Minkowski  four- 
space.  For  a  definition  see  (ll)  p.  32. 

2  (11),  p.  54. 
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is  the  "total  angular  momentum"  density.  is  the 

"canonical  energy -momentum  tensor"  : 


/ 

ik  = 


N 

w 


(F 


iFlA 


l 


(1-7) 


The  tensors 


l 

Jk^m  ~  Tki*m 


-  TLxi 


’k^m  3a 


A. 


d«C 


i,k 
N 


m 


3a 


A 


m,k 


(F,  „A  -  F.  k0) 
\  iso  ™  km  l' 


(1-8) 


may  be  interpreted  as  the  "orbital  angular  momentum"  density 

and  "intrinsic"  or  "spin  angular  momentum"  density  respectively. 

L,  n  and  S,  „  do  not  separately  satisfy  equations  of  the 
kiJm  k#  m 

form  (1-5),  hence  the  orbital  angular  momentum 


and  the  spin 


are  not  conserved  quantities,  that  is  their  values  depend  on 
the  choice  of  2.  Nevertheless,  if  2  is  taken  to  be  the 
hypersurface  =  const.,  then  the  spin  is 
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(1-9) 


This  is  identical  in  form  to  the  corresponding  expression 
for  a  free  Maxwell  field  (ll). 


The  orbital  angular  momentum  may  be  re -defined  to 
satisfy  a  conservation  law  by  using  the  symmetrised  in 
place  of  the  canonical  energy -momentum  tensor  in  (1-8) : 


Jk^m  "  Tki*m  Tkmx£ 

Since  =  T^k  ,  we  have 
definition  obscures  the  intrinsic  or 
angular  momentum  density,  since  from 


J.  „  .  =  0.  Such  a 
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that  is  the  tensor  S 


cancels  identically. 
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For  static  solutions  of  the  field  equations,  (1-4) 


reduces  to 


tip,p  0 


from  which  it  follows  that 


x  T,  .  Q  d3X  =  0 
7  iP,P 


(1-10) 


Using  (^Tlp)jp  =  67pT1(3  +  x7TlpjP  and  the  four-form  of 
Gauss'  theorem,  we  find  that  static  solutions  of  (i-l) 


satisfy  the  integral  relations: 


f  Tl7  d3x  =  j  ^  Tlpdsp  . 


(1-11) 


In  particular  static  solutions  for  which  X  T.^  -**  0  as 


X. 


7 


oo  satisfy 


T±7  d3x  =  0  . 


The  Lagrangian 
L 


=  JilV  Aljk,  tf)  d3x 


of  the  system  is  related  to  its  rest  mass: 


The  "rest  mass"  M  is  simply  the  total  energy  W  of 
the  system,  evaluated  in  the  rest  frame. 
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This 


M  =  -  I  T42|  d^x 


relation,  as  is  shown  in  Appendix  B 


n  /  ax 

X7  3t 


M  =  -L  -  77TT  /  -rr  *  ET  d^X  - 


N 

XtF 


,  can  be  expressed  as 
.  (1-12) 


. 


II.  PLANE  WAVE  SOLUTIONS 


Although  the  existence  of  solutions  of  the  field 
equations  ( I -1 )  corresponding  to  plane  waves  is  assured  by 

the  fact  that  Maxwell’s  free  field  equations  (with  the  fixed 

o 

gauge  Ak  =  -1)  are  a  special  case  of  (i-l),  the  existence 
of  non -Maxwellian  plane  wave  solutions  is  also  indicated. 

In  order  to  reduce  the  field  equations  to  a  more  tractable 
form  we  consider  solutions  having  the  particular  form 

A(r,t)  =  A(r,t)eT, 

.2  2.2 

<J>  =  y  A  s 

where  e  is  a  unit  vector  fixed  in  space  and  y  is  a  scalar 
constant.  For  such  a  case,  the  field  equations  (in  vectorial 
form)  reduce  to 

2 

e(AA  -  -  V(e-VA  +  =  [l  +  (l  -  72)A2]Ae 

St 

and  (H-1) 

y  AA  +  e*V-||  =  y [l  +  (l  -  y2)A2]A  . 

For  plane  wave  solutions,  we  assume  that 

A(r,t)  =  A(a) 

where  o  =  lc*r  -  cut,  Tc  and  being  real  vector  and  scalar 

constants  respectively.  Equations  (Il-l)  then  reduce  to 
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[e(k2-^2)  +  (a>y-k  )]  ^-§  =  [A  +  (l-72)A3]  e 

da 

2  (H-2) 

(7k2— Oik  )  =  [A  +  (1-72)A3]  7 

e  da 

where  k  =  |  1c  |  and  k  =  fce  .  The  solutions  of  (II -2) 

6 

must  also  satisfy  the  generalised  Lorentz  gauge 
3  „  [(1  +  A  2)A. ]  =  0  ,  which  takes  the  form 

OX^  K  1 

(ke-7“)  §5  [A  +  U  -  72)A3]  =  0  . 

For  nontrivial  solutions ,  this  implies  that 

k 

e 

7  =  nr  * 

whereupon  the  equations  (II -2)  reduce  to  the  single  equation 

2  k  2 

(k2-^2)  ^-4  =  A  +  (1 - ^rr)  A3  .  (II -3) 

da^  a) 


The  solutions  of  (LI -3)  may  be  expressed  in  terms  of  Jacobi 
elliptic  functions.1  If  k2  >  k02  >  o?  the  general  solution 
of  (II -3)  may  be  written  as 


A  = 


4  l+L2 


sd 


•L 


-  Tb  K  = 


7 


1+U 


where  A  and  7  are  arbitrary  real  constants,  while 

1  Definitions  and  properties  of  Jacobi  elliptic  functions 

may  be  found  in  (l4),  p.  145* 
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2  =  2.  (A2  +  -^-j)  and  L2  =  1  +  -g  g-— 

(k  /^)2-l  k2-co2  V(k2-®2), 


Correspondingly, 


-  .  ke  «  cn(Acr  s/l+L2  -  n,  K) 

E  =  dA  (“e  -  —  K)  — £- - J-Z-? - -  £ 

dn  ( Aa  y  1+L  -  r] ,  K ) 


and 

s="ixe  . 

The  relative  directions  of  the  various  vectors  is  shown  in 
figure  2.1. 


Figure  2 . 1 
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The  invariants1  of  the  field  are: 

S’-TT  =  0 

2 

and  E2  -  IT2  =  E2  (l  -  %  [l  -  cos2  (Ic, ST)]}  ,  which  can  be 

(XT  ^ 

>  or  <  0  depending  on  the  value  of  cos  (k^ET),  the  cosine 
of  the  angle  between  the  vectors  Ic  and  S’;  E  =  I  S’  | 

The  corresponding  charge  and  current  densities  are : 


P 


U2(72-l) 


1] 


and  7  =7pe  , 

where  =  (T#  ip) 

The  energy -momentum  tensor  takes  the  form 
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with  coordinates  chosen  as  indicated  in  figure  2.1. 
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Reference  (12),  p.  63* 
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l 

W  =  ^  [E2  +  (^)2  -  (M)2  +  (1+72)a2  +  |(1-72)(1+372)A4  +  K] 
p  =  [E21c  -  (f-lc)f  +  [l  +  (1-72)A2)  7A2e  ] 

We  note  that  unlike  for  Maxwellian  plane  waves  IT,  I?  and  k 
do  not  form  an  orthogonal  triad.  The  electric  field  always 
has  a  component  in  the  longitudinal  direction  (parallel  to 
k) .  In  fact  if  k2  =  k  2,  £,  e  and  £  become  collinear  and 
H  =  0  so  that  a  purely  electric  longitudinal  plane  wave 
results.  Also  unlike  Maxwellian  plane  waves,  there  exist 
non-null  charge  and  current  densities.  Furthermore,  the 
current  density  "j"  is  generally  not  parallel  to  the  direction 
of  propagation  of  the  wave. 

When  a)2  >  k2  >  k  2  the  situation  is  analogous 

J~1“  0# 

to  the  previous  one  except  that  the  solutions  of  ( II —3 ) 
now  take  the  form 

A  =  [L  cn  (j^cj  -  K) 

p  2  (jD2-k2  1  - - - 

where  \i  =  2(Mj)  — p - p  an<^  =  2  +  .  2  /  2  2N 

-k  2  K  (or -IT) 

e 


1 


The  constant  N  which  appears  in  the  definition  (1-2) 

of  £  is  taken  to  be  +1  as  for  Maxwell's  free  field  equations. 
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Therefore 

£  =  (kelc  -  to2?)  sn  (|a  -TV„K)  driver  -  rj,K) 
and  the  invariant 

2 

E2-H2  =  E2(l  -  K?)  +  )2  is  always  greater 

CCT 

than  zero.  For  the  special  case  7=0  (hence  k0  =  0) 

IiT,  it  and  Ic  form  an  orthogonal  triad  of  vectors;  nevertheless 
there  still  exists  a  non-null  current  density  T  parallel 
to  S'  ,  that  is  in  a  transverse  direction. 

2  2  2 

Finally,  the  case  where  k  >  ^  >  kg  does  not 

correspond  to  any  conceivable  physical  situation j  since  the 
solutions  of  (II-3)  are  then  reciprocals  of  Jacobi  elliptic 
functions  and  are  therefore  periodically  singular. 

The  above  solutions  are  not  generalisations  of 
Maxwell’s  free  field  solutions  in  the  sense  that  it  is  not 
possible  to  choose  the  parameter  7  in  such  a  way  that 
(II -1)  reduce  to  Maxwell’s  free  field  equations.  For  such  a 
possibility  to  arise  it  would  be  necessary  to  introduce  some 
parameter  into  the  expression  1  +  A^.2  in  such  a  way  that 
q  +  j\.  2  0  for  some  particular  value  (or  values)  of  this 

parameter.  This  is  impossible  if  the  universal  constant 
m  /  0.  Thus  the  nonlinearity  cannot  in  any  sense  be  considered 
as  a  perturbation  of  the  Maxwellian  case. 

1  Some  other  nonphysical  solutions  are  discussed  briefly  in 

Appendix  C. 
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The  field  equations  (i-l)  do  not  have  solutions 

which  would  correspond  to  spherical  waves,  since  this  would 

require  that  the  "argument  function"  a  take  on  the  form 

a  =  f (r)  +  g(t).  The  field  equations,  in  spherical  polar 

coordinates'*'  contain  terms  of  the  form 
2  2 

d  -t—  =  4^-  *  $§-  *  — i  ,  which  (for  nontrivial  solutions) 
orot  dr  dt  d(J2 

implies  that  the  product  ^  ^  must  be  a  function  of 

a  only  and  this  is  possible  only  if  o  is  of  the  form 
cr  =  kr  -  ^t .  With  such  a  o  however  the  equations  do  not 
reduce  to  ordinary  equations  in  a. 


The  general  case  of  the  field  equations,  for 
It  /  o,  was  investigated  at  some  length,  but  without  any  success. 

For  the  purpose  of  investigating  possible  particle - 

like  (i.e.,  localised)  solutions  it  is  best  to  use  spheroidal 

2 

coordinates  of  one  sort  or  another.  However,  even  in  the 
simplest  spheroidal  coordinates  -  the  spherical  polar  coordinates 
the  field  equations  constitute  a  very  complicated  system  of 
four  nonlinear  partial  differential  equations.  This  system 
of  equations  is  exhibited  in  Appendix  D.  The  difficulty  is 
compounded  by  the  fact  that  the  generalised  Lorentz  gauge 

1  The  field  equations  in  spherical  polar  coordinates  are 
given  in  Appendix  D.  Throughout  this  thesis  (r,  0,  cp)  will 
be  used  to  represent  spherical  polar  coordinates. 

2  The  more  familiar  ellipsoidal  and  spheroidal  coordinates 
are  discussed  in  Reference  (l 6),  Chapter  5* 


'6  Si  j  Ui"  ■  tv  \  :  .■  it  i  \6&  :•;  oct  IrKiji  ■  :aoo  Mx/ow  tfo'jtrtw 


■  v>  . 


•*::>  n  viar-T;  .f  M  9-  :  si  iupsrc 


. 


V 

!  sr.cit«j;oa  Xs"i>  •  >n  *<//■)  do.tdw  *  4-^/  ■  ~ 

■ 


. 


' 

r .  *»  v  ;$ 

' 

1  •  <4 

- 

i  ■ . .  a  •'..■•  •  .■:•  i  .v.  .  •-  .  ; 


/ 


. 


i  r  /t-  •awoHi  *  .‘terfsti.e  10  3rioa  g>rtc>  lb  sstarrJtbiooo 


j 

■ 

1  1  ' 


■ 


"  '  :  '  ■  ■  , 

. 


.. 


$ r 


n 


19 


which  is  inherent  in  the  theory1  provides  very  little  leeway  for 
simplifying  this  system  of  equations.  Attempts  were  made  at 
reducing  this  system  of  equations  to  a  more  tractable  form 
by  considering  various  special  cases ,  where  at  most  two  com- 
ponents  of  A±  are  non-zero  and  either  or  =  0. 

In  all  cases  where  reduction  to  reasonably  tractable  equations 
was  possible,  the  results  were  found  to  be  either  trivial 
or  unphysical. 

1  See  Chapter  I,  p.  5* 
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III.  STATIC,  PARTICLELIKE  SOLUTIONS  WITH  H  =  Q. 


The  field  equations  (i-l)  take  on  a  simple  form 

if  the  magnetic  field  strength  H  =  0.  For  the  corresponding 

dA± 


static  case,  when 


cTt" 


=  0,  they  reduce  to 


and 


(1  +  A2  -  <*>2)  A  =  0 


A<t>  =  (1  +  f?  -  <t>2)4> 


(III-D 


(III -2) 


The  corresponding  Lagrangian  (from  (1-2)  )  is 


L  =  - 


N 

Wfr 


[_  1  (v<!>)^  +  |(X2-<J>2)  +  ^(t2-<P2f+  K]  d3x. 

(Ill -3) 


Equation  (HI-l)  implies  that  1=0  or/and  1  +  1^-^  =  0. 

If  1  +  s2-*2  =  0  everywhere  then  equation  (III -2)  just 
becomes  Laplace’s  equation  which,  as  is  well  known,  possesses 
no  particlelike  solutions  (except  the  trivial  solution 
<P  =  const.)*  II  1=0  everywhere,  then  (III -2)  implies  that 
the  scalar  potential  $  must  satisfy  the  equation 


A4>  =  <J>  -  $  . 


(Ill -4) 


1 .  The  constants  N  and  K  are  chosen  so  that  the  energy 
is  finite  and  positive,  and  in  agreement  with  Maxwell’s 
theory  where  1\  -  Q  +1  =  0.  See  Ref.  (7)  und  Cnapters 

VI  and  VII . 
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This  particular  situation  has  been  considered  by  Dr.  Schiff  (7), 
and  the  results  are  summarised  in  Chapter  V.  A  discrete, 
countable  set  of  particlelike  solutions  does  exist  in  this 
case*  As  is  clear  from  (ill -4),  they  are  asymptotic  to  zero 
as  ^ — -  Y ( 0 , cp) 1 *  for  large  r,  and  therefore  correspond  to 
particles  of  zero  charge  (since  ;  Ab  d3x  =  0) . 


The  most  general  implication  of  equation  (III-2) 
is  that  the  condition  1  +  1^  -  ^  -  0  hold  in  some  region 
V  of  space,  while  the  other  condition  1=0  apply  elsewhere. 
On  the  boundary  S  between  these  regions  both  conditions 
must  apply:  1=0  and  b  =  -  1.  Furthermore,  £  =  -  Vb 
must  exist  and  be  continuous  everywhere  for  a  particlelike 
solution,  hence  in  particular  on  S.  The  equation  (III -2) 
then  becomes 


Ab  =  b  -  b3  in  V,  (HI-5) 

Ab  =  0  outside  V,  (III-6) 


with  b  =  -  1  on  the  boundary  surface  S.  If  V  is  a 
closed,  simply  connected  region  of  space  then,  choosing  the 
origin  of  coordinates  inside  V,  the  general  well  behaved 
solution  of  (III-6)  may  be  written  as 


1  Y(0,cp)  is  some  linear  combination  of  the  spherical 

harmonics . 
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*  =  t  1  -  a  +  2  2  -jj2L-  Y™  (0,q>)  ,  (III -7) 

m=-i  r 

a;  b^m  being  constants  and  Y™  the  usual  spherical  harmonics. 
The  boundary  surface  S  is  represented  by  r  =  R(d,cp),  where 


00  &  b  o 

2  2  Y1^  =  a  .  (Ill -8) 

4=0  m=  " fr+1  ^ 


The  vector  potential  A,  outside  V,  is  given  by 

=  4>2-l;  its  direction  must  be  radial  if  H  is  to  be  zero 

for  any  <J>  satisfying  Laplace’s  equation.  Inside  V ,  $  must 

be  a  well  behaved  solution  of  (ill -5)  with  ^  =  -  1  and 

50  =  _  ”  |  ^+l)b.em  ~m,„  <  „  1 

^  i=o  m=-X[H(S,cp)  1 


2+5  Y^(0,(p)  on  S 


The  physical  interpretation  of  such  solutions ,  if 
they  exist.,  is  determined  by  the  values  of  the  constants 
bgm.  Thus,  if  all  bim  =  0  save  bQ0  =  b,  then  the  solution 
clearly  corresponds  to  a  spherically  symmetric  charge  particle 
of  charge  b(— ).  2  The  boundary  surface  in  this  case  is  a 

s 

1  Note  that  the  continuity  of  -g-n-  and  ,  which  are  required 
if  Vb  is  to  be  continuous  on  S,  is  assured. 

2  Since  for  every  solution  <t>  of  (III-5)*  Is  also  a 

solution,  it  is  clear  that  such  particles  occur  in 
oppositely  charged  pairs. 
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sphere  of  radius  R  =  -  *  inside  which  the  spherically  symmetri- 
cal1  charge  density  is  non-zero.  Similarly *  if  all  b^m  =  0 
except  bQ0  =  b  and  b1Q  =  p*  the  solution  corresponds  to  a 
charged  particle  with  charge  and  electric  dipole  moment 
proportional  to  b  and  p  respectively .  The  boundary  surface 
in  this  case  is  the  spheroid 


=  R(n)  =  1 5  f1  +  ./  1 


/ 


+ 


4ap 


cos 

~2 


0 


with 


<  1. 


The  charge  densities  corresponding  to  such  charged 
particlelike  solutions  are  given  by 


P 


0 


in  V 

elsewhere . 


(Ill -9) 


Generally  speaking*  the  expression  (III-7)  represents  the 

electric  multipole  expansion  of  the  particle’s  localised 

2 

charge  distribution  (III-9). 


1  It  is  shown  in  Appendix  E  that  well  behaved  solutions  of 
(III -5)  which  are  spherically  symmetric  on  S  are  such 
throughout . 


This  does  not  mean  that  particlelike  solutions  corresponding 
to  any  values  of  b«  are  possible.  For  example*  the  scheme 
does  not  admit  a  neutral  particle  with  only . a  pure . electric 
dipole  moment:  all  b»  =  0  save  b-,Q  =  p^  since  this  implies 
that  the  boundary  surface  S*  given  by  ar  =  p  cos  0*  is 
not  real . 
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It  is  possible  to  consider  more  complex  subdivisions 
of  space  into  alternating  regions  of  linearity  (where 
1  +  ^  -  p2  =  o  and  p  =  0)  and  nonlinearity  (where 

/ 

A  =  0  and  p  ^  0) .  We  consider  here  only  the  case  of  spherical 
symmetry;  assume  a  subdivision  of  space  into  spherical 
regions  by  concentric  spheres  of  radii  0  <  R-^  R2  ^  ^ 

where  <D(R±)  =  ±  1,  for  all  i  =  b  2,  ....  n.  The  innermost 

region  (where  0  r  'C  R^)  is  necessarily  a  nonlinear  region, 
since  otherwise  the  only  well  behaved  solutions  are  the 
trivial  solutions  0  =  -  1. 

If  n  =  1,  we  have  the  situation  discussed  above, 
n  =  2  implies  that 

AP  =  0  -  0^  for  0  <  r  <  R1  and  R2  <  r  <  » 

=  0  for  R^  <  r  <  R2 

with  0  =  -  1  on  r  =  R-^  and  0  =  ^:1  oft  ft  =  Rg* 

Well  behaved  solutions  of  this  system  correspond  to  neutral 
particles . 1  Such  solutions  are  discussed  in  Chapter  VI. 

1  As  is  pointed  out  in  Chapter  IV,  the  only  well  behaved 
solutions  of- (III -4)  in  R2  <  r  <  «  are  those  asymptotic 
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Similarly,  when  n  =  3*  well  behaved  solutions  of 
the  system  of  equations 


Ab  =  b  -  b^ 

if 

0  <  r  <  Ri 

and 

r2 

< 

r  *  R3 

=  0 

if 

FL  <  r  <  R2 

and 

R3 

< 

r  <  oo  , 

with  b  =  -  1  when  r  =  R-^,R^  an<^  b  =  +  1  when  r  -  Rg, 

correspond  to  charged  particles,  since  in  the  outermost 

region,  R^  ^  r,  the  well  behaved  solutions  are  of  the  form 

a  +  Such  solutions  are  discussed  in  Chapter  VII. 
r 

In  general,  for  any  n  0,  well  behaved  solutions 
(when  they  exist)  correspond  to  neutral  particles  when  the 
outermost  region,  R^  <  r  <  00  ,  is  "nonlinear"  and  to 
charged  particles  when  it  is  "linear". 
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IV.  EXISTENCE  AND  UNIQUENESS  OF  SOLUTIONS 

•3 

OF  A0  =  0  -  0 


It  is  evident  from  the  discussion  of  the  previous 
chapter  that  the  existence  of  neutral  and  charged  particlelike 
solutions  (with  ft  =  0)  is  closely  tied  in  with  the  existence 
of  well  behaved  solutions  of  the  nonlinear*  elliptic  partial 
differential  equation 

o  1. 

Ab  =  b  -  (IV-1) 

Phase  space  analysis  of  the  equation  for  the  case  of  spherical 

symmetry  (15)  suggests  the  existence  of  a  countable  infinity 

of  discrete  global  solutions  which  are  asymptotic  to  zero  as 
-r 

£ _  *  together  with  a  continuum  of  solutions  asymptotic  to 

r  _ 

+  i  as  sin(y/2  r  +  ft)  ^  p  being  a  constant.  The  existence 

of  the  discrete  set  of  solutions  asymptotic  to  zero  is 

suggested  also  by  the  fact  that  equation  (IV-1)  may  be 

2 

written  in  the  "Schrttdinger"  form: 

[_  A  -  01 2]  0  =  -0  , 

1.  The  equation  possesses  three  trivial  solutions  0  =  0*  -  1. 

2.  This  was  pointed  out  by  Dr.  Schiff*  (7)* 
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which  may  be  interpreted  as  the  wave  equation  for  the 
eigenfunctions  b  of  a  particle  with  fixed  "energy"  -1, 
moving  under  the  influence  of  an  attractive  potential 
V(r)  =  -b2.  This  point  of  view  suggests  the  existence  of 
non-spherical  eigensolutions  in  addition  to  those  of 


spherical 

symmetry. 

The  fact  that  well  behaved  solutions  of  (IV-1) 

do  in  fact  exist  is  easily  demonstrated:  For  the  sake  of 

simplicity  we  consider  solutions  of  (IV-1)  which  are  axially 
.  ^  0  _  « 

symmetric  (^r  =  0). 

The  substitutions  f  (r,  p.)  =  rb ( r,  p) 
and  p  =  cos  0  reduce  (IV-1)  to 


2  3 

if  +  K  P[f]  -  f  +  ^  =  0  ,  (IV-2) 

dr  r  r 

where 

p[f]  s  |r  |r  ]  • 

Assuming  that  solutions  of  (IV-2)  are  well  behaved. 

let 

f(r,n)  =  2  cc  (m.)  ri+1  (IV-3) 
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Then 


2  00  00  n 

=  2  a0([i){£+l  )S  r^"1  =  2  ( S+3)  (S+2)  a »  2  (p)  r^+1  , 

dr  i/=0  1  i=  -1  • 


^  P[f] 
r 


f3 

r 


00 


i 

2  P 
=0 

Cai(l 

00 

00 

00 

2 

2 

2 

m=0 

n=0 

00 

i 

m 

2 

2 

2 

x=o 

m=0 

n=0 

J-l 


00 


=  2  P  [a^2(M,)]  r 


S+l 


S>=-2 


and 


Og(l-0  am(^)  an(p)  r 


i+m+n+1 


a„  (p)  a  (p)  a  (p)  r 
i-rrr^'  m-nv^'  nv^' 


i+1 


Substituting  these  expressions  into  (IV-2)  and  equating  the 
coefficients  of  the  various  powers  of  r  y  we  find  that 


a  (n)  =  c  ,  an  arbitrary  constant 
0  o 


Pla^p)  3  +  2*1  ax(p)  = 


=  0 


P  [a,  p((i,)]  +  (^+3)  (i+2)a.  2(p)  o 

(IV-4) 

JL  m 

“  alM  -mlo  nio  aH- m^Vn^  anM  ’ 


for  i=0,  1,  2,  ... 

This  is  a  set  of  linear,  but  inhomogeneous  Legendre  equations 
and  the  existence  of  well  behaved  solutions  (for  I  p  I  <  l)  is 
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easily  established  by  induction.1  The  general  well  behaved 
solution  of  (IV-4)  may  be  written  as 

^(p-)  =  +  Co3  °1 3  ck-2^  3 

where  Pk  is  the  Legendre  polynomial  of  order  k,  ck  an 
arbitrary  constant  and  yk  a  polynomial  of  degree  k-2  in 
p,.2  The  general  solution  for  ak(p)  thus  contains  k-1 
arbitrary  constants. 

It  is  shown  in  Appendix  G  that 

Og(-|±)  =  (-1)^  Og( |l) 

hence,  as  a  corollary,  the  following  theorem  applies  to 
the  solutions : 

f  (r,  -p)  =  -  f (r, p)  if  and  only  if  a^(p)  =  0  for 
f  odd  1 

all  i  \  &.* 

L even  J 


1  This  is  done  in  Appendix  F. 

2  It  is  a  special  solution  of  the  complete  equation  for 

ak(p).  A  definition  of  yk  is  given  in  Appendix  P. 
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It  now  remains  to  consider  the  radius  of  convergence 
of  the  indicated  series  solution.  We  consider  the  case  of 

-p 

spherical  symmetry*  for  which  =  0  or  equivalently 

c  =  0  for  all  k  >  0.  The  coefficients  a  are  now 
k  K 

constants  (which  depend  on  c  )  and  as  such  are  trivially 
even  functions  of  d,  so  that  =  0  for  all  odd  4.  For 
even  4  *  we  obtain  from  (IV-4) 


(4+ 3) (4+2)  o ti+2 


4 

2 


m=0 


m 

2 

n=0 


0Ln  CL  CL 
4-m  m-n  n 


(IV-5) 


Using  the  ratio  test*  it  is  easily  seen  by  induction  that 


the  series  converges  uniformly  to  an  analytic  function  in 
the  region  0  <  r  <  R(cq)  for  every  finite  cQ  *  although  R 


decreases  with  increasing 


In  the  region 


^  <  r  <  R’*  where  R’  is  arbitrarily  greater  than  R*  the 


equation 


d2f 

^7 


f  + 


r 


=  0 


(iv-6) 


now  satisfies  a  Lipshitz  condition*  hence  Picard’s  method  may 

3 

be  used  to  establish  the  existence  of  a  unique  solution. 


1  This  is  established  in  Appendix  H. 

2  Some  explicit  values  of  f(r)  as  computed  from  the  series 
have  been  worked  out  for  various  aQ.  They  are  listed  in 
Appendix  Jf. 

3  (14),  p.  83,  189. 
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Thus  for  the  case  of  spherical  symmetry,  (IV-1)  has  a 

unique  global  solution  for  each  value  of  4>(r=0).  All 

d.45 

these  solutions  have  the  property  that  =  0  at  r  =  0 
and  4>(-r)  =  4>(r).  Furthermore  phase  space  analysis  of  the 
equation  (IV-6)  (Reference  15)  implies  that  all  these 
solutions  are  finite  everywhere. 

For  the  more  general  case  of  axial  symmetry,  it 
is  similarly  seen  that  the  series  solution  (IV-3)  converges 
in  a  finite  region  about  the  origin.  The  general  series 
solution  contains  an  infinite  set  of  arbitrary  constants 
c,  ,  for  k  =  0,  1,  2,  ...  It  is  not  obvious  however  that 

K 

these  constants  may  be  so  chosen  that  the  series  converge 
globally  to  nontrivial  functions. 

The  series  solutions  exhibited  here  are  useful 
only  for  determining  the  behaviour  of  the  solution  in  the 
vicinity  of  the  origin.  It  is  impossible  to  infer  from  them 
how  global  solutions  of  (IV-1)  behave  for  large  values  of  r. 
For  the  case  of  spherical  symmetry  this  information  may  be 
obtained  with  the  aid  of  phase  space  analysis,  while  for  more 
general  cases  it  is  not  even  obvious  that  global  solutions 
exist . 

The  problem  of  determining  charged  particlelike 
solutions  takes  on  the  form  of  a  Cauchy  problem  for  the 
nonlinear  equation,  in  that  well  behaved  solutions  of  (IV-1) 
are  sought  in  a  closed  region  V  of  space,  with  2  =  -  1  and 
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given  everywhere  on  the  bounding  surface  S.  Since  the 
equation  (IV-1)  is  elliptic  one  might  surmise  that  the 
problem  is  poorly  posed  (overspecified) .  In  fact  the 
boundary  value  (Dirichlet)  problem  is  also  poorly  posed 
(underspecified)  in  that  if,  for  example,  S  is  taken 
(arbitrarily)  to  be  the  sphere  of  radius  R  =  2.3,  then  there 
exists  a  discrete  set  of  solutions  (in  all  probability  an 
infinity  of  them,  corresponding  to 

r1 2-.|l  =  0,  .123,  2-706,  6.279,  12.138,  20.161,  29.591,  etc., 
1 

on  r  =  2.3 . 

As  already  remarked  the  situation  here  resembles 
most  closely  (at  least  for  the  spherically  symm.  case) 
the  eigenvalue  problem  for  the  Schrtfdinger  equation,  rather 
than  either  the  Cauchy  or  Dirichlet  problems.  Unlike  the 
Schr&dinger  problem,  however,  the  eigenvalue  does  not  appear 
explicitly  in  the  equation  but  is  rather  contained  in  the 
boundary  conditions . 

For  nonspherical  charged  particles,  however,  the 
boundary  conditions  correspond  to  Cauchy  data,  hence  the 

2 

existence  of  the  "interior"  solution  (reV)  is  rather  obscure. 

1  The  first  few  (numerically  obtained)  solutions  are  plotted 
in  fig.  4.1. 

2  As  is  the  concept  of  corresponding  to  an  eigenvalue 

since  is  no  longer  a  constant  on  S. 
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Figure  4 . 1 
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V.  APPROXIMATION  METHODS  FOR  PARTICLELIKE  SOLUTIONS 


Nonlinear  field  equations  of  the  type  considered 
here  generally  cannot  be  solved  explicitly  by  means  of  the 
usual  techniques;  various  approximation  methods  must  be  used. 
One  of  the  simplest  and  most  successful  of  these,  especially 
for  obtaining  approximations  to  eigensolutions ,  is  the 
variational  method: 


Consider  field  equations  of  the  form 


A4>  =  Fl(0) 


(V-l) 


I  d  P  1 

where  F  ($)  =  ^  and  F  is  some  simple  function  of  F. 
Since  well  behaved  solutions  of  (V-l)  also  extremise  the 
Lagrangian 


L(4>) 


r 


(V4>)2  +  F(4>)  ] 


(V-2) 


this  suggests  that  approximations  to  the  solutions  of  (V-l) 
can  be  obtained  variationally  by  extremising  L  with  respect 
to  suitable  comparison  functions. 


I  S 

1  For  example,  if  F  (<D)  =  b  -  in  some  localised  region  V 

of  space  and  zero  elsewhere,  then  (V-l)  is  identical  with 

(III -5). 
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If  0  is  a  solution  of  (V-l)  and 

o 

some  suitable  variational  trial  function  (u 
but  small  perturbation),  consider  the  nature 


f  =  4)  +  u 

being  an  arbitrary 

of  the  bound 


L(f)  to  L(<t>) : 


Kt)  = 


=  L('t’)  +  k-  ([W-Vu  +  u  Fl(<l>)]d3x 


r 


+  w 


+  hr 


[^(Vu)2  +  |  u2  F"(<>)]d3x  (V-3) 

S  j"  XT  F^b'1>)d3x  . 


Using  V*  (uV4>)  =  Vu*Vd  +  u  A4>  and  Gauss’  theorem: 


r 


[  V0  •  Vu  +  u  FT(0)]d^x  =  0  uV0  •ds'  + 


[ -A4>  +  Ft (0) ]u  d^x 


Since  u  is  arbitrary  this  implies  that  L  is  extremised 
if  and  only  if  d  is  a  solution  of  (V-l).1 


1  Note  that  for  solutions  which  do  not  vanish  at  infinity 
the  indicated  surface  integral  does  not  in  general 
vanish  and  then  the  variational  principle  which  is 
equivalent  to  the  field  equations  is 

5L  _  |  5<t>  VcD-ds*  =  0  rather  than  5L  =  0.  For  example, 

this  is  the  case  for  charged  particlelike  solutions 
discussed  in  Chapter  VII. 

Throughout  this  thesis  the  terms  "variational  trial  function" 
and  "variational  comparison  function"  will  be  used  inter¬ 
changeably. 
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If  u  is  sufficiently  small,  the  nature  of  the 
variational  bound  is  determined  by  the  second  order  term 


62L  =  ^  j  [|(Vu)2  +  |  u2  f"  (♦)  ]chx.  (V-3a) 

M 

Since  F  (<t>)  is  in  general  indefinite,  so  presumably,  is 
the  variational  bound  L(^)  to  L(4>).  Suppose,  however,  that 
the  variational  comparison  functions  f  are  restricted  to 
those  which  satisfy  the  integral  relation 

f  1* 

[(WO2  +  6  F(^)]d3x  =  0  ,  (V-4) 


which  relation  is  also  satisfied  by  the  solutions  of  (V-l). 
Replacing  f  by  ^  +  u  in  (V-4),  we  get,  to  second  order  in  u 


tj  u2P  ($)d3x 


[—V^  *  Vu  +  ^(Vu)2  +  u  F,(4))]d3x  . 


Thus,  to  second  order  in  u  (V-3)  becomes 


L(^)  =  L(«>)  +  ^ 


r 


J  Vb-Vu  d3x  + 


±-(Vu)2d3x.  (V-5) 


Since  $  is  a  solution  of  (V-l)  it  extremises  L  hence  the 
first  order  term  in  (V-5)  vanishes,  while  the  second  order 


1.  This  is  equivalent  to  requiring  that  the  Lagrangian 

be  an  extremum  with  respect  to  a  radial  scale  parameter 

of  an  arbitrary  comparison  function  ip,  (p.  40 )  .  • 
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term  is  now  positive  definite.  Thus,  if  the  variational 
comparison  function  ip  is  chosen  to  satisfy  (V-4)  then  L(f) 
provides  an  upper  hound  to  L/^)  provided  \  ip  -  $  I  is  suffi¬ 
ciently  small.  This  differs  significantly  from  the  usual  linear 
eigenvalue  problem  in  that  no  other  conditions  (such  as  the 
subsidiary  orthogonality  condition  in  the  linear  problem)  are 
necessary  to  ensure  upper  bounds  to  excited  states. 

The  condition  (V-4)  imposed  here  on  the  comparison 
function  ip  is  possibly  not  the  only  one  which,  by  itself, 
can  assure  that  L (if/)  give  an  upper  bound  to  L(4>)  for 
|  f  -  b  |  sufficiently  small.  It  is  clear  that  any  condition 
satisfied  by  the  exact  solutions  of  (V-l)  may  be  used,  pro¬ 
vided  it  may  be  utilised  to  eliminate  the  generally  indefinite 

O  Q  p  * 

term  /  F"  (b)u  dDx  from  5  L,  and  reduce  the  latter  to 

a  positive  definite  form. 


Note  that  in  general  the  integral  relation  (V-6)  below 
(which  is  equivalent  to  extremising  L  with  respect  to  the 

amplitude  of  ip)  does  not  by  itself  assure  that  L(^)  is  an 

2. 

upper  bound  to  L($). 


1.  It  may  happen  that  F" (b)  is  a  positive  definite  quantity,  in 
which  case  L (ip)  will  give  an  upper  bound  to  L(<t>)  for  any  ip 
such  that  \p  -  ^  is  small.  This  is  not  the  case  for 

F  (0)  =  <t>  -  <J>3. 

2.  This  is  because  when  (V-6)  is  solved  (to  second  order  in  u) 
for  /F"X<J))u1 2  d^x  one  of  the  terms  which  arises,  namely 
-1/2  /  u^FMI( b)d3x,  is  also  indefinite  in  general.  However, 
as  pointed  out  by  Dr.  Schiff,  for  the  case  FT(b)  =  $  -  $ 

it  is  positive  definite. 
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In  using  the  variational  principle  to  obtain 
approximations  to  the  solutions  of  (V-l),  the  variational 
Lagrangian  may  be  extremised  with  respect  to  variable 
amplitude  and  scale  parameters  without  otherwise  specifying 
the  nature  of  the  comparison  function. 

Let  ^(x,y,z)  =  A  f(ax,  Py,  yz,  6i) 

where  A,  (a,  (3,7) ,  (i  =  1,  2,  ...)  are  amplitude,  scale 
and  other  parameters  while  f  is  some  suitable  and  integrable1 
comparison  function.  For  L(f)  to  be  an  extremum  with  respect 
to  A  we  require  that 

Iff-  w  W-  Is  W  +H]d3x  =  0 

Since  (Vji)  =  V(||)  =  \  W  and  |f  =  F>  (f)  ||  =  F'W 

The  above  condition  reduces  to 

/** 

[  (v^/)2  +  ^F’(^)]  d^x  =  0.  (v-6) 


1  "Integrability"  here  implies  that  all  indicated  integrals 

involving  f  (hence  f)  exist. 
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The  nature  of  this  extremum  is  determined  by  the  value  of 

.2  n2  . 

— £  at  the  "extremal"  value  of  A.  Noting  that  — K-  =  0, 


_1 _ 

4ttA2 


[(W)2  +  if/2  F"(f)  ]  d3x  , 


hence  at  the  extremum  (where  (V-6)  holds): 


d2L 


$)A 


4tta‘ 


[if/2  F"  (if/)  -if/  F1  (if/)  ]  d3x  . 


ext . 


Thus,  whether  the  extremum  is  a  maximum  or  a  minimum  (or 
neither)  depends  on  the  explicit  form  of  F (if/).  In  the 

O 

particular  case  F1  (if/)  =  if/  -  if/D 

= - o  \ d3x  <  0, 

4ttA  3 

ext . 


so  that  L  is  a  maximum  with  respect  to  A. 

To  extremise  L  with  respect  to  scale  parameters, 
consider  the  transformation  x  =  —  ,  whereupon  the  Lagrangian 
may  be  written  (with  p  =  y  =  l) 
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L 


1 


ri  f^(x',y,z)x2  1 

l2  1  dx  >  +  2 


(||) 2  +  |(||)2  +  F(^)]dx  dy  dz 


r 


(H>2 


+ 


2a 


s^2  +i 


y 


2a 


dy  dz 


Therefore 


r 


Ba 


1 

w 


1 


Mi2  _  .  _i_  i^t)2  - 1 


2a 


2  v  ay ' 


2  a 


2  v  dz 


■^r  F(^) ]dx Tdy  dz  =  0 


a 


implies  (after  performing  the  inverse  transformation  xT  =  ax): 


r 


r 


(H)2  A  - 


j 


[  (V*)2  +  2  F(f)  ]  d3x 


Analogous  results  hold  for  y  and  z;  hence 


r 


r 


2  d3x  = 


'ax' 


r 


(g)2  A  - 


r 


d3X  = 


(k^t)2  +  P(V')  ]d3x 
(V-7) 


if  L  is  to  he  an  extremum  with  respect  to  scale  parameters. 
Relation  (V-4)  follows  immediately  from  (V-7 ) . 


As  already  mentioned  the  integral  relation  (V-4) 
is  also  obtainable  by  extremising  L  with  respect  to  a 
radial  scale  parameter.  This  particular  extremum  is 
actually  a  maximum  for  any  F(^),  which  implies  that  integrable 
(hence  particlelike)  solutions  of  (V-l)  are  unstable  with 
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respect  to  a  uniform  radial  "stretching"  of  the  particle.1 

Finally  we  note  that  the  integral  relations  (V-4), 
(V-6)  and  (V-7)  must  be  satisfied,  in  particular,  by  integrable 
solutions  of  (V-l). 

In  using  the  variational  method  the  usual  procedure 
is  to  choose  a  comparison  function  which  is  explicit  in  its 
dependence  on  the  independent  variables  but  contains  also 
variable  parameters.  For  example  a  comparison  function  for 
a  possible  odd-parity  eigensolution  of  (III-4)  having  axial 
symmetry,  might  be 

b  =  A[r  e ~ar  P  1(p)  +  b  r^  e"^r  P^^p.)],  (V-8) 

where  and  P^  are  Legendre  polynomials,  while 

A,  a,  f3,  7  and  b  are  variational  parameters.  Substitution 

of  (V-8)  into  the  Lagrangian  then  yields 

L  =  L(A,  a,  p,  7,  b) 

which  can  then  be  extremised  with  respect  to  the  parameters, 
"optimal"  values  of  which  are  thus  obtained. 

1  This  was  pointed  out  by  R.H.  Hobart  (Proc.  Phys .  Soc./ 

London/  82,  201  (1963))  for  spherically  symmetric  solutions 
and  by  Derrick  (l8)  for  general  integrable  solutions 
of  (V-l). 
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It  is  possible,  however,  to  choose  comparison 
functions  for  which  the  dependence  on  one  (or  more)  of  the 
independent  variables  is  not  specified.  To  illustrate  the 
point,  consider  a  variational  comparison  function  of  the  form 

<t  =  e’ar  Y(n),  (V-9) 

where  a  is  a  variable  parameter  and  Y(p)  is  some  function 
to  be  determined.  The  Lagrangian  then  becomes 

L(Yja)  =J  oC  (Y,Y*,M*;a)  dp.  . 

The  resulting  Euler -Lag range  equation  will  yield  an  ordinary 
differential  equation  which  may  be  solved  (numerically,  if 
necessary)  for  the  optimal  form  of  Y(d)  and  hence  L.  It 
should  be  noted  that  even  this  "optimised”  L  will  (in 
general)  still  contain  variable  parameters  (a  in  the  example 
above)  and  will  therefore  have  to  be  extremised  with  respect 
to  them. 

Approximations  to  particlelike  solutions  of  (V-l) 
can  also  be  obtained  numerically.  In  the  case  of  spherical 
symmetry  (when  (V-l)  reduces  to  an  ordinary  differential 
equation)  a  number  of  methods  exist  for  obtaining  numerical 
solutions,1  of  which  the  Runge-Kutta  and  finite  difference 


1  See,  for  example.  Fox  (19) 9  Chapter  I. 


}  o;-  ■ 


,r  '' 

-  V 


' 

*  '  •  •  ■  XT'"  •. 


ri  < 


:  .  .....  .  I  .  •-•.o'-  i  .  J  '  •  "H  t'-:  : 

(P.-V) 


J 


-m)  Y  "i: 


. 


.  >oc(  no  U  '  '  *n;  i  -aJ  -sriT  .i:  $c:  s.L-.b  '-d  ot 


(.'  •  t.-j ■  ,'fetY.)  Oo 


-  -  i . 


.  1  a  j  v  '  •:  i  •  •  •  '• : ;  t" 


. 

' 

r  -  ' 


. 

■  ’ '  Is  =•'  '-VJ 


[s^-j  .  '  ..  n  s  t  - 


: 

>  * 


. 

■ 


ii>3.  '  o  *  1£,0  t  '!  •  f.  i  >;■•  a  'it  o?  i.  •:  :■ 


„  .  • 

■ 


■  . ... >  ;•  .  ’  : ..  ■  o.  ;  S-X!0 i cUlCG  .  ' 


— 


. 


i  ’■ 


A 


...  - 


43 


methods  are  best  known.  The  latter  method  has  been  used 
to  obtain  all  numerical  solutions  quoted  throughout  this 
thesis  (except  where  stated  otherwise).1 2  The  finite 

difference  analogue  to  (IV-6)  is 

f(rH-5)  +  r(r-8)  -  2f(r)  _  f(r)  +  fRuL  0)  (V-10)  , 

r 


where  5  is  a  small  increment  in  r.  In  replacing  (IV-6) 

4 

2  dr 

by  (V-10)  terms  of  the  order  of  5  — jr  and  higher  are 

dr 

neglected.  Approximate  particlelike  solutions  of  (IV-6) 

may  then  be  obtained  by  iteratively  solving  (V-10) .  For 

example  for  neutral  particles,  the  asymptotic  form  of  the 

solution  (for  r  greater  than  some  arbitrary,  large  R)  is 
-r 

A  - —  .  Values  of  the  solution  for  r  less  than  R  can 

r 

then  be  computed  from  (V-10)  for  various  A.  Actual 
approximate  solutions  of  (IV-6)  then  correspond  to  those 
values  of  the  "asymptotic  amplitude"  A  for  which  the  behaviour 
of  the  numerically  computed  f(r)  near  the  origin,  is  of  the 

p 

required  form.  The  iteration  can  also  be  performed  in 


1  All  ordinary  differential  equations  considered  here  are 

reducible  to  the  form  ,2  which,  according  to 

S-S  =  f(x,y) 
dx 

Fox  (19,  p.  4)  is  most  easily  and  accurately  solved  by  the 
finite  difference  technique. 

2  The  form  that  the  particlelike  solutions  of  (IV-6)  must 
have  in  the  vicinity  of  the  origin  was  discussed  in 
Chapter  IV  (  (IV-3)  and  (IV-5)  ) • 
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the  reverse  order  (starting  near  r  =  0)  as  was  done  (for 
neutral  particles)  by  Tesbima  (8). 

A  similar  procedure  may  (hopefully)  be  used  to 
obtain  numerical  approximations  to  non-spherical  particle¬ 
like  solutions,  if  such  exist.  The  field  equation  (III-4), 
in  spherical  polar  coordinates,  is 


1  8 
"IT  c>r 


/  2  , 

(r  + 


r  sin  0 


H - p - p— 

r  sin  0 


8  /  .  a 
(sin  6  se^ 


^ 


(V-ll) 


=  0 


8cp 


The  behaviour  of  all  particlelike  solutions  in  the  vicinity 
of  the  origin  was  determined  in  the  previous  chapter  (  (IV-3) 
and  (IV-4)  ). 


For  charged  particles  the  behaviour  of  solutions 
in  the  "linear"  region  (where  -  1  =  0)  is  given  by 

(III -7).  For  neutral  particlelike  solutions,  the  asymptotic 
form  for  large  r  is  (from  (V-ll),  with  |  <t>  |«  l) 


d 


00  l 
2  2 
jG=0  m=-j0 


te  (e'cp) 


where  the  Y™  are  the  usual  spherical  harmonics,  while 

kg(r)  =  -i^  h^1'^  (ir) ,  h^(x)  being  spherical  Hankel  functions 

of  the  first  kind  (  (17*  P*  539)  )•  Since 
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-r 

k.(r)  -*  - —  as  r  ->  oo  (for  all  &)  then,  for  large  r, 

Jo  r 

solutions  of  (V-ll)  which  are  asymptotic  to  zero  are  of  the 

-r 

form  4>  =  ^ —  Y(0,<p),  where  Y(0,cp)  is  some  linear  combina¬ 

tion  of  the  spherical  harmonics.  In  any  case  the  behaviour 
of  particlelike  solutions  for  large  r  is  readily  determined. 
Thus,  if  odd-parity  neutral  particlelike  solutions  of  (V-ll), 
having  axial  symmetry  are  sought  then,  for  r  greater  than 
some  arbitrary  but  large  R  the  asymptotic  form 


$  =  2  A^k^(r)  P2i+l^cos  may  be  used*  For  r  ^  R  and 

£=0 

0  <  G  <  tt,  we  replace  (V-ll)  by  its  finite  difference  analogue 

c3  <t>  .  * 

First  however,  since  the  coefficient  of  in  (V-ll)  is 
singular  at  0  =  0  ,  it  is  useful  to  make  the  transformation 
[i  =  cos  0  (and  f  =  r4>),  whereupon  (V-ll)  becomes 


+  ^ 


b  r  /_  2  s  Sf  i 

b;:  [(lHi  ) 


f1 * 3 

f  +  ^ 

r 


0. 


(V-12) 


Its  finite  difference  analogue  is 


1  In  analogy  to  the  transformation  f  =  r4>,  there  exists 
a  further  transformation  which  eliminates 

f  =  - 

v^sin  Q 

the  first  partial  with  respect  to  0  (or  |i)  from  (V-12)  . 
However,  since  this  transformation  brings  in  singular 
coefficients  (at  0  =  0)  into  the  resulting  equation  it  is 
not  desirable  to  carry  it  through  (at  least  if  the 
equation  is  to  be  solved  numerically). 
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f  ( r+5  ,  p.)  +  f(r-5,p.)  -  2f(r,p.)  +  1-pT  j  f  (r,p.+A)+f  (r,p.-A)  -2f (r,p 


2m- 
r 


:  (r,„;  A)-f  (r,„  a)|  _  f  (r^)  +  f3(^)  =  o  , 


(v-13) 


with  terms  of  the  order  of  6 


2  d^f 


A2  d^f 


o^f 


TT  '  T  7T’  2  7T 

or  r  op.  r  op. 


and 


smaller  being  neglected  (6  and  A  are  small  increments  in 
r  and  p.  respectively)  .  Thus  (V-13)  might  be  used  to  compute 
values  of  f  for  r  <  R  and  |  p.  |  <  1 .  A  solution  would 
then  correspond  to  those  values  of  the  constants  hr  for 
which  the  behaviour  of  f  near  the  origin  is  of  the  form 
required.  Clearly,  this  procedure  for  nonspherical  states  is 
much  more  cumbersome  since  the  asymptotic  form  of  the  solution 
now  is  indeterminate  to  within  an  arbitrary  analytic  function 


00 


P(p.)  =  2  AiPiId). 

1=0 


At  p.  =  -  1  "forward"  and  "backward"  difference  approxi^ 
mations  (rather  than  central)  must  be  used  for  df  and  o^f 

dp. 

since  f  is  undefined  for  |  p.  |  >  1. 


dp/ 


2  The  numerical  methods  described  here  were  used  successfully 
(as  a  trial)  to  reproduce  the  known  fs  and  p  state 
eigensolutions  of  the  Schrttdinger  equation  with  the 
potentials  V(r)  =  -i.  and  -e-^'. 


I 

/'  0  V'- 


a:  ■ 


. 


v  7 


1 


, 


. 


. 

. 


I  ■  ' 

■  > 

;  .  -  ' '  '*■  • 

r  *  (*)'* 

- -mm** 

■r 


.  I 

. 


. ;  <  |  i,  |-  tot  .Tear  or  ai 


- 


.  ,  -t  \  -  •  1  i.'  '  i'C£>3  '£,*•  '  .  CO  .  -■  ,  ■  ■  V  ;:-n  '«  - 

bna  J.*  »  (i)V 


VI.  NEUTRAL  PARTICLES 


Neutral  particlelike  solutions  of  the  static  field 

equations  (III -4)  were  obtained  variationally  by  Betts , 

Scbiff  and  Strickfadden  (10)  and  numerically  (for  the  case 

of  spherical  symmetry)  by  Teshima  (8).  In  agreement  with  the 

implication  of  phase  space  analysis  (15) ,  numerical  solution 

of  (III-4)  yields  a  discrete  set  {0,  3  of  neutral  particle- 

1  -r 

like  solutions  asymptotic  to  zero  as  A^  — —  .  The  first 
four  are  plotted  in  figures  6.1  to  6.4. 

As  shown  by  Betts  et  al  (10),  these  eigensolutions 
satisfy  the  "orthogonality"  relation 

f  0.^  b.  d^x  =  /0.^0.  d^x  (VI-1) 

J  i  J  J  J  i 

Since  1  =  0,  the  spin1 2  of  the  particles  represented  by  these 
solutions  is  zero,  as  expected.  Also,  since  it  =  0,  the 
Poynting  vector  fT  =  ^  S'  x  IT  is  zero,  hence  so  is  the 

o 

angular  momentum  J^m  (or  L ^m) .  The  Lagrangian  of  the  system 

is,  according  to  (III-3) , 

L  =  -  -^r  /[-^(V0)  +  ^  0  -  jj-  d^x  (VI -2) 

1  Recall  that  the  "spin  density"  is  proportional  to  1  x  t. 

2  These  were  defined  in  Chapter  I. 
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Fourth  Neutral  Particlelike  Eigensolution 
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Note  that  the  constants  N  and  K  which  appear  in  (III-3) 
have  been  chosen  to  be  -1  and  0  respectively,,  to  assure  a 
finite, positive  mass  spectrum.  The  sign  of  the  Lagrangian 
for  neutral  particles  is  thus  opposite  to  that  of  Maxwell’s 
theory,  hence  also  opposite  to  that  for  the  charged  particle 
system.  This,  as  pointed  out  by  Dr.  Schiff  (j) ,  implies 
that  charged  and  neutral  particles  are  distinctly  separate 
systems  within  this  scheme . 

The  masses  of  the  particles,  according  to  (1-12) 
(and  using  the  integral  relations  (V-4)  and  ( V-6)  ),  may 
be  written  as 

M  =  ^  /  ** 1 2d3x  =  /  <2  d3x  =  U-jr  /  (V4>)2  d3x. 

These  values  for  the  first  few  neutral  particlelike  solutions 
1 

are  : 


(n) 

A 

Vo) 

M 

(0) 

-  2.713 

*4.34 

1.503 

(1) 

■t  17.11 

+14.11 

9.46 

(2) 

i  83.87 

-29.15 

28.71 

(3) 

-375.0 

749.43 

63.97 

.  .  J n  ,  .  .  lim  r , 

1  (n)  is  the  number  of  nodes  of  the  solution.  A  =  ’  re  0 . 

2  3* 

A,  0(0)  and  M  are  in  units  of  —  ,  and  mc g 


g 


3 


g 


respectively . 
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It  is  reasonable  to  assume  that  there  exist  also 
solutions  of  (III-4)  which  are  not  spherically  symmetric.1 2 
Variational  approximations  to  a  possible  non-spherical 
eigensolution  of  the  form  <t>(r,p)  =  -<J> ( r,  -p ) ,  that  is  an 
axially  symmetric  solution  with  a  nodal  plane,  were  obtained 
by  Betts,  Schiff  and  Strickfadden  (10)  using  the  method  of 
parameter  variation.  Their  results,  for  the  lowest  (in 
energy)  odd  parity  state,  are  reproduced  in  table  VI -A. 

They  seem  to  suggest  the  existence  of  an  odd  parity,  eigen¬ 
solution  corresponding  to  a  neutral  particle  of  rest  mass 

<  o  mc^tf 

~  3  •  9o  — —  * 

S 

The  variation  of  function  method,  described  in 
Chapter  V,  was  used  by  this  author  to  obtain  approximations 
to  possible  odd  parity  eigenstates.  This  method  requires 
quite  a  bit  more  numerical  computation  than  parameter  variation. 
Extremisation  of  the  Lagrangian  with  respect  to  the  undetermined 
functions  which  are  contained  in  the  chosen  trial  function 
yields  a  set  of  ordinary  differential  equations  which  have  to 


1  Some  arguments  in  favour  of  such  an  assumption  are 
discussed  by  Dr.  Schiff  (7)  and  in  Chapter  IV. 

2  This  method  is  described  in  Chapter  V. 
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Table  VI -A.  Variational  Approximations  to  Lowest  Odd-Parity 

Axially  Symmetric  State. 


Comparison  function 

Optimal  values  of 

parameters 

Mass 

A 

a 

P 

d 

0 

b 

1.  A  z  e'ar 

32 

7T 

/3 

5.47 

-ar^1 

2  .  A  z  e 

3.03 

2.7  6 

1.8 

4.72 

—  2?  * 

3.  A  z'  e 

32 

3 

/5 

V 

T 

3 

4.41 

4  .  A  fg 

.833 

4.14 

5.  A[p1(n')r’e'r' 

+P3(n'  )r'?enir’ 

+b  P3(M.')r5e'CTr] 

6.48 

2.08 

2.18 

1.8844 

2.807 

.6939 

3.96 

In  the  table : 

rT  s  (xf,  yr,  zT)  =  (ax3ay3p>z) ,  rT 

1  rT| 

j  dT 

_  z  T 
r T 

:  Legendre  polynomial  of  degree  n. 

me 

Mass  is  in  units  of  — g—  . 

g 

The  function  f  was  obtained  by  Strickfadden  (9)  as  a  solution  of 

s 

Ab  -  b  +  b^  =  0  using  Greens  function  methods.  0T  is  the 
comparison  function  no.  1  of  the  table.  The  results  for  com¬ 
parison  function  no.  5  (previously  unreported)  are  due  to 


Dr.  Schiff. 


V  \> 

v 

• 

- s'  • 

-  *V  . 

-  X- 


I . 


>  -  ’  -  .  I  V  •M-'.-l  .  •.  S 


r 


— 


♦  s.ijs.i  p.  \l  *3.  'A 

■ 


—  .  «I  -L. 


■  * 

■ 

■  II  !■*>>■ 


■■l08 . 


- 


■  :'b8.  r. 


■  j 


r  - —  "  i1”’ 


Fv 

„  ■ 

Jr.3 

|  g . 


■> 


—— 


- 


A 


•vr 


v 

50,  e 

2' 


££S. 


8*i.S 


•10- 


,  -  jo  —  ■< 

>S 


■ 

e  ’  s  'A  .8 


r  '  «* 


1  ' 


.  f 


t 


t  r:t  > 


^  -  *lo  acting  rt±  at  aasM 


S  ■ 


' 


•Mi;  at  Tt>  .  afccrUom  nol  ron*Jl  anso^tO  antaxj,s  0  * 

i 

o:.t  ojjJb  9*i£  (bs^noqornu  ^laaoivs'iqt)  $  ,  on  rrot$onj/i  nc  . 


•  T.  L.X:  •••  ' 


55 


be  solved  numerically.  The  procedure  does  have  the  advantage 

that  the  set  of  trial  functions  so  obtained  satisfy  identically 

some  of  the  integral  relations  which  are  satisfied  by  the 

q  1 

exact  solutions  of  -  r.  *  In  particular  they 

satisfy  the  "orthogonality"  relation  (VI-l),  suggesting 
that  the  various  "excited"  trial  eigensolutions  actually 
furnish  approximations  to  corresponding  "excited"  nonspherical 
states  of  the  field  equation,  if  such  exist. 

The  following  types  of  trial  functions  were 
considered : 

4>A  =  ft(r)  P1(n) 

A  m  -ar  ,r/  \ 

<$>B  =  r  e  Y  (p) 

$c  =  R(r)  P  1{[i)  +  Q(r)  (ijl) 

=  (r1 2  +  b  r^)e  ar  Y(p) 

^E=  R(r)  P1(p)  +  Q(r)P3(p)  +  S(r)P5(n) 

2 

where  P  (p,)  are  Legendre  polynomials,  a,  b  and  m  -  parameters 
and  R,  Q,  S  and  Y  functions  for  variation.  Details  of  the 
calculations  together  with  graphs  of  the  (numerically  obtained) 

1  This  is  shown  for  each  of  the  trial  functions  considered 
in  Appendices  K,  L,  M  and  0. 

2  m  takes  on  integral  values  only. 
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optimal  forms  of  the  functions  R,  Q,  S  and  Y  are  given  in 
Appendices  K,  L,  M,  N  and  0.  The  results  for  the  first  few 
odd-parity  eigenstates  are  listed  in  Table  VTvB. 


units  of 


In  the  table  An , 

—  ,  while  M  is 
g 


Bn,  Cn  and  Y 
in  units  of 


(n)w=w 

mc^f) 

2 

g 


are  in 
The  integer 


n  is  used  to  label  the  eigensolution  in  order  of  increasing 
M.  (Except  in  the  case  of  ^or  n  >  2,  n-1  is  equal 

to  the  number  of  nodes  of  the  functions  R^n\  S^, 

and  Y^n^ .  ) 


The  nodeless  optimised  variational  trial  functions 
presented  in  table  VI -B  seem  to  approximate  the  lowest  two- 
particle  state1  rather  than  a  possible  single  particle  state 
of  odd  parity,  whose  existence  is  suggested  by  the  results 
presented  in  table  VT-A  (especially  the  trial  function  5) • 


1  vis.  <t>  2=  VI  ^  I)  -  ^0(!  ^  ^  where  <1>0  is  the  lowest 

spherically  symmetric  one-pjrticle  state  (p.  52)  and  |  T? 
Recall  that  M(4>0)  =  1.503  hence  M(*12)  =  3-006  S—s?1 

(Reference  (id)). 
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Table  VI -B.  Variational  Approximations  to  Odd  Parity 

Neutral  Particlelike  Solutions. 


A.  cj^n)=  R(n)  (r)Px(d)  [R^n^  (r)  ^  Ank1(r)  as  r  -> 


n  AM 

n 


1 

±  -10.08 

4.87 

2 

±--73.64 

l6 . 60 

3 

±394.82 

37.88 

B. 

-©■ 

“a 

m 

=  r 

e-ar  Y<n>  (m.) 

n 

(|-L=1 )  m 

a 

M 

1 

±  36.66  2 

2.19 

3.76 

2 

±  19.59  2 

1.32 

31.33 

3 

±10.70  2 

0.874 

109.8 

C. 

<t,(n) 

c 

=  R^n 

) (r)P1(d)+Q^n^ (r)P3(d) 

[ R^ n  ^ ( 

r)->Anki(r) 

q(d)  (r)-*Bnk3(r)  as 


n 

1 

2 

3 

4 


±  10.15  1  .8618  3*63 
±  80.00  ±  25.75  10.72 
±  70.46  ?  27.34  26.75 
±276.15  *  63.46  38.92 


D. 


E. 


=  (r2  +  br3)  e-“r  Y<n)(b) 


n 

1 


Y^n  ^ (|±=l) 


+ 


30.47 


a 

1.76 


b 

-0.228 


3 


<t> 


(n)  =  R(n)  (r)pi([JL)  +  ( r)  P^  (|i )  +  S  ^  (r)  P^M-) 


[R 


E 
(n) 


Anki(r)*  Q 


(n) 


B  k_(r)  and  S 
n 


(n) 


cn*5(r) 


n 


M 


11.30  ±  1.4137  ±.048532 


r  -►  00  ] 


M 

.73 


as  r  -►  00  ] 


1 


3.31 
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Attempts  were  made  to  obtain  odd-parity  eigensolutions 
of  (III-4)  using  the  numerical  technique  described  in  Chapter  V. 
The  difference  equation  (V-13)  corresponding  to  (ill -4) 
was  solved  by  iteration  assuming  that  the  eigensolutions 
have  the  asymptotic  forms 

a)  4)1  =  A  k1(r)P1(p)  1 2 3 

and 

b)  <t>2  =  A1k1(r)P1(q)  +  A^ ( r)  P3 (q)  . 

For  case  a)  this  was  done  using  A  =  .1,  .2,  .  ..,  50.  No 

meaningful  results  were  obtained.  It  was  found  that  as 

r  -►  0,  the  "solutions"  oscillate  violently  for  all  values 

of  A  considered;  the  amplitude  of  these  "oscillations" 

2 

increases  with  decreasing  r. 

A  number  of  reasons  for  such  behaviour  come  to  mind: 

1.  No  well  behaved,  global^  odd-parity  solutions  of  equation 
(III-4)  exist,  which  are  asymptotic  to  zero.  The  results 

3 

of  variational  calculations,  however,  suggest  the  contrary. 

1  P  are  Legendre  Polynomials  while  k^  are  spherical  Hankel 
functions  with  imaginary  argument  (see  p.  44). 

2  A  sample  result,  for  A=10  is  listed  in  Appendix  P,  together 

with  plots  of  (r,  p)  for  some  (fixed)  values  of  r. 

3  See  also  the  arguments  given  in  Chapter  IV,  and  (7). 
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2.  The  well  behaved  solutions  are  "obscured".  In  the 
numerical  solution,  by  the  possible  existence  of  solutions 
which  are  unstable  and/or  divergent  near  the  origin. 

Since  no  such  solutions  exist  in  the  case  of  spherical 
symmetry,  it  would  be  curious  indeed  if  such  appeared 

in  the  odd-parity  case. 

1. 

3.  No  solutions  exist  with  the  asymptotic  form  4^ 

4.  The  iterative  method  of  solving  the  finite  difference 
approximation  to  (III-4)  which  is  used  here,  is  not  an 
efficient  and  accurate  method  for  this  type  of  equation. 
Some  comments  on  the  method  are  given  in  Chapter  V. 

As  noted  there,  the  method  has  been  used  successfully 
for  certain  linear  odd  parity  eigenvalue  problems. 


If  the  third  possibility  is  taken  to  be  the  most 

probable,  then  a  more  realistic  asymptotic  form  that  b1 

is  necessary.  This  may  be  obtained  from  the  variational 

comparison  functions  4  of  table  VI -A  and  of  table  VT-B, 

which  are  of  the  form  <t>2  indicated  on  p.  58  .  In  each  case 

2 

only  a  small  region  of  the  Ax,  A^  plane  was  investigated  in 


1.  It  will  be  recalled  that  the  general  asymptotic  form  of 

the  solutions  is  00 

b  =  Z  A.k„(r)P*(|i)  . 
i=0  Ji  ib  1 


2.  It  is  clear  that  the  amount  of  machine 
investigate  exhaustively  even  the  case 


time  necessary  to 
of  two  asymptotic 


parameters  A*  (say  A^  and 
such  an  undertaking. 
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the  vicinity  of  those  values  for  which  corresponds  to 

the  indicated  trial  functions.  The  results  obtained  were  in 
all  cases  qualitatively  similar  to  those  obtained  with  the 
asymptotic  form  $1  .  Finally,  the  variational  trial  function 
5  of  table  VI -A  was  used  to  approximate  an  odd  parity  solution 
in  the  asymptotic  region.  In  this  case  also  the  numerical 
results  obtained  could  not  be  interpreted  as  approximating 
a  well  behaved  solution. 

As  noted  in  Chapter  III,  neutral  particlelike 
solutions  are  possible  when  space  is  subdivided  into  alter¬ 
nating  regions  of  linearity,  where 


and 


A4>  =  0 


-  1 


1 


and  nonlinearity,  where 


and 


A4>  = 


1=0. 


The  linear  regions  are  bounded  by  concentric  spheres  of 
radii  Rn  and  Rfi+1  (0  <  Rn  ^  Rn+1  ^  00 with  =  -  1  on  R, 


1.  Note  that  |  l'\  is  imaginary  since  I  |  <  1 
region . 

2.  We  consider  here  only  spherically  symmetric 


in  the  linear 

solutions . 
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1 

and  0  =  +  1  on  Rn+1  ^  so  that  every  linear  region  must 
contain  a  zero  (node)  of  the  solution.  Thus  a  solution 
possessing  3  nodes,  say,  can  be  one  with  0,  1,  2  or  3 
regions  of  linearity.  There  are,  then,  a  total  of  eight 
possible  distinct  neutral  particlelike  solutions  having 
3  nodes.  In  general,  if  (N)  denotes  the  class  of  solutions 
having  N  zeros,  then  the  greatest  number  of  distinct  neutral 
particlelike  solutions  which  can  belong  to  this  class  is 

,  .  S  N(N-l)  ...  (N-k) 

1  +  - kl 


The  mass  corresponding  to  a  given  particlelike 

2 

solution  is,  from  (III-3) 


M  = 


d3X  +  2  i-  /  i(V$)2d3x, 


where  2  and  2  denote  summation  over  all  linear  regions 
J  nil 

Vg,  and  all  nonlinear  regions  ^  respectively. 

1  In  addition  ^  must  be  continuous  on  Rn  and  Rn+1 • 

2  Recall  that  N=-l,  K=0  and  M=-L  for  neutral  particlelike 
solutions . 
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This  expression  can  be  reduced.,  with  the  aid  of  the  integral 
relation  (VIII-2),  to 


1 

T6rr 


) 


implying  that  the  linear  regions  contribute  nothing  to  the 
rest  mass  of  the  particle. 


In  the  table  below  we  list  the  (numerically  obtained) 


results  for  some  compound  neutral  particlelike  solutions : 


1 


Type 

A 

4(0) 

1  ;L) 

±  17.01 

+  14.15 

2;L,NL) 

±  83.85 

±  29.15 

2 jNL,L) 

±  80.90 

±  29.42 

2;L,L) 

±  80.89 

±  29.42 

3;L,NL,NL) 

±375.0 

+  49.43 

3 ;NL, L,NL) 

±374.3 

•f  49.47 

3jL,L,NL) 

±374.3 

+  49.47 

3 ;NL,NL,L) 

±342.2 

+  50.18 

3;L,nl,l) 

±342.2 

•f  50.18 

3;  NLj  L, L) 

±341.7 

■f  50.21 

3;L,L,L) 

±341.7 

■f  50.21 

In  the  table 

A,  0(0)  and  M  are 

in  units 

M 


9 .46 
28.71 
28.58 

28.58 

63.97 

63.96 

63.96 

63.44 

63.44 

63.43 

63.43 


tie  me' 


g'  S 


and 


3 

mc respectively.  (2;NL,L)  refers  to  a  solution  with 
S 

2  zeros,  with  the  region  around  the  first  (where  |  4>  |  <  l) 
being  nonlinear  and  around  the  second  -  a  linear  one. 
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These  results  differ  very  little  from  those  for  the  corres- 

1 

ponding  solutions  listed  on  page  52.  This  is  because  the 
regions  of  linearity  (at  least  for  the  first  few  solutions 
considered  here)  are  near  the  origin  and/or  small  in 
extension  and  |  <t>  |  <  1  there.  In  fact  the  only  perceptible 
difference  (to  3  significant  figures)  in  corresponding  mass 
values  occurs  where  there  is  a  linear  region  at  some  distance 
from  the  origin^  (for  example  (2;NL,L)  or  (3;NL,NL,L)  ). 

A  sample  curve,  (2jNL,L),  showing  the  linear  and  nonlinear 
regions,  is  plotted  in  figure  6.5* 


1  For  this  reason  they  are  not  plotted,  since  on  a  scale 
like  that  used  in  figures  6.1  to  6.4  they  are  indistin¬ 
guishable  . 

2  The  region  near  the  origin  contributes  little  to  the  mass 
because  of  the  factor  r1 2  in  the  integrand  of 
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VI I .  CHARGED  PARTI CLES 


In  this  section  we  consider  static  charged  particle¬ 
like  solutions  of  the  field  equations  in  the  absence  of  the 
magnetic  field.1  The  simplest  such  solutions  are 

Ai  =  (ffj  10)  when  0  <  r  <  ^ 

and 

=  (±  Jb2-1  fj  i(a  +  p)  )  when  r  >  -  , 

2  A 

where  a  and  b  are  constants  ,  r  is  a  unit  vector  in  the 
radial  direction  and  0  a  solution  of: 

^  (r2  ^)  =  f  -  4>3  when  0  <  r  <  |  ,  (VII -l) 

r 

with  q 

2  J 
.  ,  d0  a  ,  b 

0  =  -1  and  when  r  -  -  . 

1  See  Chapter  III,  p.  21. 

2  We  shall  use  both  a  and  a  =  ±1  -  a  throughout  this 
chapter . 

3  Since  for  every  solution  0,  -0  is  also  a  solution  it  is 

sufficient  to  consider  0  =  -1  on  r  =  b/a  instead  of 

b  =  ±1 . 
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The  Lagrangian  associated  with  each  solution  is, 
from  (III-3), 

L  =  /  [|  (|f)2  +  ^2  -  ^  ^  -  if]r2dr  +  \  cb.  (VII-2) 

0 

The  constant  N,  which  appears  in  (III-3),  is  taken  to  be 
unity  so  that  the  energy -momentum  tensor,  in  the  region  where 
1  +  Ai2  =  0,  will  be  the  same  as  in  the  Maxwell  theory. 

K  is  taken  to  be  ^  so  that  the  Lagrangian  and  hence  the 
energy,  which  is  linearly  related  to  it1,  is  finite.  The 
corresponding  mass  is  then 

M  =  -L  -  ab. 

2 

The  spin  and  orbital  angular  momentum  are  zero 
as  for  neutral  particles. 

As  shown  in  Chapter  IV,  there  is  a  solution  of 
(VII -1)  for  each  value  of  <D(0)  (the  constants  a  and  b,  of 
course,  change  with  b(0)  )  .3  These  solutions  may  be  grouped 

1  Equation  (1-12) . 

2  These  were  defined  in  Chapter  I. 

3  b(o)  is  the  amplitude  of  b  at  the  origin. 
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together  according  to  the  number  of  nodes  they  possess  on 
the  lines  0=0  and  0  =  -1 .  We  shall  use  the  notation 
(M,  N)  to  denote  the  group  of  solutions  of  (VII-l)  having 

1 

M  nodes  on  the  line  0=0  and  N  nodes  on  the  line  0  =  -1 . 

There  exists  *  then,  a  continuum  of  well  behaved 

solutions  of  (VTI'-l)  within  each  group  (M,  N)  .  They  may  be 

characterised  by  a  single  scalar  parameter  such  as  0(0)  or 

2 

b  (the  charge  parameter). 

A  discrete  set  of  particlelike  solutions  can  be 
determined,  for  example,  by  imposing  the  requirement  that 
the  charge  parameters  be  all  the  same  (in  absolute  value)  or 
have  small  integral  ratios.  It  is  preferable,  however,  to 
use  a  general  principle  to  determine  a  discrete  set  of 
particlelike  solutions,  rather  than  an  artificially  imposed 
one  like  that  cited  above.  For  this  purpose  the  principle  of 
least  action  may  be  invoked:  Those  solutions  for  which  L 
is  a  minimum  with  respect  to  the  parameter  (or  parameters) 
characterising  the  continuum  (that  is,  the  most  stable  solutions) 
will  be  considered  as  representing  charged  particles. 

1  Some  such  solutions  are  plotted  in  figures  7 .2  -  7*5- 

2  The  parameter  a(or  a)  is  not  convenient  for  characterising 
the  solutions,  as  is  shown  below. 
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Thus,  the  discrete  set  of  spherically  symmetric 
charged  particlelike  solutions  is  characterised  by  the 
discrete  set  1^]  ,  for  which  L(b)  is  a  local  minimum, 
and  the  mass  M±  =  -L(b  )  -  b^fb^  is  positive. 

If  b(b)  and  b(b  +  5b)  are  two  "neighbouring" 

solutions  of  (VTI-l)  (6b  is  assumed  to  be  small),  then  the 

change  in  L  corresponding  to  5b  =  b(b  +  5b)  -  ^.(b)  is, 

as  pointed  out  in  Chapter  V,  p.35* 

6L  =  irr  60  Vb'ds”  =  -  b  6a  , 

since  b  =  a  +  —  for  r  >  —  .  is  a  spherical  surface 

10  0 j 

of  radius  R  -►  00 .  In  the  limit  as  6b  -*■  0  (hence  5b  -►  0) 
this  becomes 

1 

dL  y,  da 
db  db  * 

Thus,  for  b  /  0,  L  and  a  have  concurrent  extrema. 

For  the  discrete  set  of  charged  particlelike  solutions 
of  ( VII -1 ) ,  the  action  principle 

da 
db  (0) 


1  Equivalently  dL 

db  (0 ) 
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5L  +  b  5a  =  0 

then  becomes  6L  =  0  subject  to  6a  =  0. 

This  may  be  cast  into  a  simpler  form  by  the  use 

of  a  Lagrange  multiplier;  that  is.,  the  above  statement  is 

equivalent  to  requiring  that 

6L  +  b  T  5a  =  0  (VII-3) 

where  br  is  an  undetermined  multiplier.  Since  for  exact 
solutions  of  (VII-l)  bT  =  b,  it  is  clear  that  a  variational 
calculation  based  on  (VII -3)  will  yield  approximations  to 
the  discrete  set  of  charge  particlelike  solutions,  and  the 
multiplier  bT  will  furnish  an  approximation  to  the  corres¬ 
ponding  "charge"  b.1 2  1  '  .v'ri.Mv  f  for 

L(b)  exhibits  a  minimum  for  each' group ■ of  solutions 
:(M,N);  an  example  ;of  such ..>  arburv’e 'corresponding  to  the 

group  of  solutions  (1,1)  is  plotted  in  lig'ure  7  .'1  together 

2 

with  a(b)  and  M(b)  for  this  group. 


1  If  no  subsidiary  conditions  (on  the  comparison  function)  are 
imposed  then  (VII-3)  will  yield  an  approximation  to  the  "ground 
state"  charged  particle,  that  is  to  that  solution  which  corres¬ 
ponds  to  the  smallest  of  iL(b^)]. 

2  The  behaviour  of  these  functions  for  other  groups  of  solutions 
is  qualitatively  very  similar  (except  for  the  group  (0,l)  for 
which  L  .  =0  at  b=0  and  a(b)  exhibits  no  extremum).  In  general 
M(b)  exBiBits  no  extrema. 
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L(b),  a  (b),  and  M(b)  for  the  group  (l,l) 


Fig.  7.1 
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The  first  few  spherically  symmetrical  charged  particlelike 

1 

solutions  correspond  to: 


group 

L 

b 

a 

M 

0(0) 

d,i) 

4.48 

2.76 

-2.177 

1.530 

6.98 

(1,3) 

7.46 

7.89 

-2.168 

9.65 

7.00 

(2,2) 

21.1 

8.02 

-3.234 

4.80 

-18.61 

(3,3) 

57.15 

15.70 

-4.246 

9.50 

35.68 

These  solutions 

are  plotted 

in  figures 

7.2  to 

7.5.  (Only 

that  portion  of 

solution  for 

which  r 

VI 

is  plotted  in 

each  case) . 


The  (nearly)  integral  ratios  of  the  charges  (at 
least  for  the  first  few  states)  is  worthy  of  note.  It  is 


possible  that  the  departure  from  integral  ratios  which  is 

apparent,  is  a  reflection  of  the  fact  that  these  results  are 

2 

only  approximate. 


•m  p  3^  C 

1  In  the  table  M  is  in  units  of  ,  b  in  units  of  — 

2  g * 1 2 * * *  6 

and  a  in  units  of  .  These  results  were  obtained 

§ 

numerically  by  Dr.  Schiff  using  the  Runge-Kutta  technique 
for  solving  (VII-l). 

2  Although  it  is  believed  that  all  numerical  results  given 
in  this  thesis  are  at  least  within  1%  of  the  exact  values, 
it  is  not  possible  to  determine  this  "deviation"  with 
certainty. 
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Fourth  Charged  Particlelike  Solution:  (3*3) 
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Solutions  for  which  the  difference  M  -  N  is  odd 
correspond  to  negative  mass  states,  hence  cannot  represent 
actual  particles.  Since  for  such  solutions  the  vector 
potential  takes  on  imaginary  values  when  |  |  <  1  it  might 

be  conjectured  that  only  solutions  with  A  real  represent 
particles . ^ 

Approximations  to  the  charged  particlelike  solutions 
may  also  be  obtained  variationally  using  the  action  principle 

5L  +  b  6a  =  0. 

To  facilitate  the  choice  of  variational  comparison  functions 
it  is  useful  to  make  the  transformation  <f>  =  X  -  1*  Equation 
(VII-l)  then  becomes 

\  (p2  |£)  +  2X  -  3X2  +  X3  =  0  ,  ^  0  <  r  <  | 

r  (vii -4) 

n  dx  a2  ,  b 

with  x=  0  and  =  -  —  when  r  =  -  . 

Correspondingly,  the  Lagrangian  is 
b/a 

L  =  f  [|  (g^)2  -  X2  +  X3  -  ^  Xb  r2dr  +  |  ab  .  (VII-5) 

0 


1  This  was  pointed  out  by  Dr.  Schiff  (7)* 
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Let  X  =  X  f(|_i*L  r^),  X,  V->  being  amplitude,  radial  scale 

and  other  parameters  for  variation,  and  f  an  arbitrary 
integrable  comparison  function  with  a  zero  in  the  region  r  >  0. 

v 

From  X  =  0  at  r  -  -  we  immediately  get 


(vn-6) 


where,  for  a  given  f,  c  is  a  function  of  r\±  only, 

defined  by  f(c,r]i)  =  0  .  Similarly,  from 
2  h 

^  when  r  =  —  we  get 

dr  b  a  ° 


v(T)±)  = 


(VII -7) 


where  the  function  v (ru  )  =  - 


df  (x,ti,  ) 


ox 


-1  x  =  c(ri1)  . 

Because  of  the  relations  ( VTI -6)  and  (VII-7 ) 
all  but  (any)  two  of  the  parameters  a,  b,  X,  p,  are 

(initially)  independent.  We  shall  take  a,  p,  T[±  to  be  the 
independent  parameters,  so  that 

L  =  L(a,M.,r]i)  , 


h(a,LL,r|p= 


ac(ii1) 

d 


^  )  v(rj1  )o(rll)  ' 
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Then  6L  -  b  5a 

-  ('5a  '  b)6a  +  5JT  6ti  +  2  5t^  5T|i  ~  0 

implies 

§  =  b  (vII-8) 

§=°  (raI-9) 

=  0  .  (VII -10) 

The  Lagrangian  (VTI-5)  now  simplifies  to 

2 

L(cx,t-,'ni)  =  G0(\)  -  ^3  t^Z|^(T]1)-A3G3('ni)+A2G2('ni)  ] 

M* 


where 

G0(ri1 )  =  o3v2  +  G(ni) 

and 

°(\)  =  f°Q  (I) 2x2  dx 

G,  (rj . )  =  /  [f(x,T).  )]k  x2  dx  for  k=2,3  and  4 

K  1  o 

Prom  (VII-8) 


V2  ■ 

=  6[-^-A2G2j  -  AG^  ^2  ^  for  A  /  Oj  (VII  ll) 

and  from  (VII-9) 

G|r2  =  A2G^  -  37G3  +  2G2 


(VII -12) 


t 


■y  . 


. 


) 


■ ,  y 


■ 


j 


•i 


i. .  ?}  foa  (ci-liv)  <ifi£3flJ3’X3£j 


% 


"i 


■ 

, 

■ 


1 


■  «•; 


■wertw 


0 


V  •' 


*  \t:  .  )  I  61 


-• 


•\r-£ 


v  -  ■  •  ••  ■•- 

Y  ~  (  *P ) 

0 


. 


sQS  ,  .  r  *{.  -  ,ES/C 


-  i|0 


, 


79 


These  relations  may  be  solved  for  A  and  \i  as  functions 
of  T) .  .  Implicitly: 


r£  1  o  n  1/2 

-  [q—  (j|-  A  -  AG^  +  Gg)  1 


where 


(2c3v2-G)G^A2  -  6(c3v2-G)GqA  +  4(c3v2-2G)Gq  =  0  . 


The  corresponding  expression  for  the  Lagrangian  is 


Results  obtained  in  this  way, 

1 

are  tabulated  below: 


f  (x,ri1) 

Numerical  inter, 
gratlon  of  (VTI-1) 

1-  x1^ 

1-2x+t]1x2 

RfL  2t]-2 

1  -2x  +r)2x 

.99463 

.635 

.99017 

1 

2+1  -  ^1 

X  +1 

.0830 

cos  (x+r]1)+ri2 

2.925 

.999878 

for  various  choices  of  f. 


L 

b 

M 

a 

4.48 

2.76 

1.53 

-2.177 

no  extremum 

found 

4.76 

1.90 

.491  - 

-1.857 

4.59 

2.27 

.930 

-2.013 

4.245 

2.25 

1.10 

-1.95 

4.77 

1.90 

.  556 

-1.841 

1  The  extremisation  of  L  with  respect  to  in  was  done 

numerically  in  all  cases.  Expressions  for  rhe  integrals 

G  ,  G^  for  each  function  are  listed  in  Appendix  q. 
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The  variational  procedure  is  not  very  efficient 
for  obtaining  approximations  to  the  charged  particlelike 
solutions.  The  choice  of  suitable  comparison  functions  is 
restricted  by  the  following: 

(1)  f(x,T^)  must  possess  a  zero  for  x  >  0, 

(2)  c  must  be  explicitly  solvable  in  terms  of 

(3)  the  integrals  Gq,  G^  must  be  expressible  in 
terms  of  elementary  functions. 

The  computations  are  rather  involved,  and  while 

the  results  for  L.  are  not  unreasonable,  the  approximate 

values  of  the  physically  meaningful  quantities  M  and  b 

are  in  poor  agreement  with  those  obtained  from  numerical 

1 

solution  of  VII -1.  It  is  clear  that  the  variational 
procedure  would  be  much  more  cumbersome  for  obtaining 
approximations  to  any  non-spherical  charged  particlelike 
solutions . 

We  consider,  finally,  "compound"  spherically 
symmetric  charged  particlelike  solutions,  that  is  solutions 


1  This  is,  of  course,  because  it  is  L  that  is  extremised 
in  the  variational  procedure  not  M  or  b.  Thus,  if  the 
variational  comparison  function  ip  -  0  +  e(<J>  being  the 
exact  solution)  then  |  L(^)  -  L(^)  |  is  of  the  order  of 

whereas  for  M  or  b  the  corresponding  deviation  is 
of  the  order  of  e  .  This  differs  from  the  neutral  particle 
case  where  L  =  -M. 
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where  the  regions  in  the  vicinity  of  their  zeros  (and  where 
|  0  |  <  l)  are  linear  regions  such  that 


and 


AO  =  0 


01 2  -  1  <  0 


there . 


The  Lagrangian  in  this  case  becomes 


L  =  s  i—  X  [|(v<t’)2+i<l>2  -  ■K*  d3x  +  z  w  /v  |(Vlt>)2d3x> 

^  *ni 


where  g  and  2  denote  summation  over  all  regions  of 

^  n  i  -L 

linearity  and  nonlinearity  respectively.  In  analogy  to 

the  neutral  particle  case  this  expression  may  be  reduced, 

with  the  aid  of  (VIII-2)2,  to 


L 


=  2 
n£ 


1 

TF 


IT 


I 


V. 


(^-l)d3: 


-  Tj-  ab. 


nZ 


Again,  only  those  solutions  for  which  =  0  and  the 

corresponding  mass  is  positive,  are  taken  to  represent 
charged  particles. 


1  Recall  (Chapter  III)  that  the  innermost  region  must  be  a 
nonlinear  one  while  the  outermost  region  (r  >  b/a)  must  be 
a  linear  one  for  charged  particles. 

2  Note  that  P(0)  =  tj02  -  ^  ^  for  charged  particles 

in  nonlinear  regions. 


. 


0  A 


j  - 


,  e.  :i 

■ 


i 


is  ,  2  xvi.-. ..  i  ritf  'i  aj30  lie  ix^q ,  XiS'io  -^n'  ©rl^ 


03  ,  (S-XIIV)  lo  fcfjB  9ri3  (13  lw 


.*»  |  -  ...  \ 


%d!  :  l,x*«  £ 


•  let.  o  .  :o  ;  i/Ios  ■  ac  "  I  no  <ni  ... 


:>....  '  •'  -  ;  ;  :<  ;  ;  ? .  ;oq  so  o  -  • 


#vB9lO  ‘Bq  b  a../1  :  • 


- 


iaorrrrsctuo  s>11  si  1.1.’  sr  ?  fnsi-'jc  tlnon 


*  ^  =  (  '  cJsn  -  o.r oM’  S  v. 


82 


A  summary  of  the  first  few  (numerically  obtained) 
solutions  is  given  below1: 


group 

L 

b 

a 

M 

0(0) 

(2,2;L) 

21.13 

8.01 

-3.234 

4.77 

-18 . 63 

(3,3;L,NL) 

57.15 

15.69 

-4.246 

9.48 

35.69 

(3,3;NL,L) 

57-09 

15. 66 

-4.238 

9.30 

35.77 

(3,3jL,L) 

57-09 

15.67 

-4.238 

9.31 

35.76 

These  "compound"  solutions  differ  very  little  from  those 
listed  on  p.  71  j  hence  they  are  not  separately  plotted. 

As  an  example  the  solution  (3j3;NL,L)  is  plotted  in  figure  7 .6. 

The  search  for  nonspherical  charged  particlelike 
solutions  was  confined  only  to  the  possible  existence  of 
charged  particles  with  a  pure  electric  dipole  moment;  that 
is,  well  behaved  solutions  of: 


1  (-3*3  jNL,L)  refers  to  a  solution  belonging  to  the  group 

(3,3)  with  the  region  about  its  first  zero  being  nonlinear 
while  that  around  its  second  zero  being  linear.  The  group 
(l,3,L)  contains  no  particlelike  solution. 


2  L  and  M 


units  of 


me 


are 

2 


in  units  of 


mc^fi 


a 


and 


is  in  units  of 


and 
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A0  =  4>  -  $3  when  r  <  R(|i) 

(VII-13) 

0  =  ti_a  +  ^+  ^  for  r  >  R(q), 

r 

where  b  and  p  are  constants  proportional  to  the  charge 
and  electric  dipole  moment  respectively.  The  surface  of 
the  charge  distribution  is  given  by 

r  =  R(m-)  -!=[*+(!  +  ^)1/Z  ] 

b 

b2 

where  \i  =  cos  0  ;  for  a  real  surface  |  p  |  <  I  7^  I 
is  required. 

Attempts  at  obtaining  well  behaved  solutions  of 
(VII -13)  numerically,  using  the  same  method  as  for  non- 
spherical  neutral  particlelike  solutions,1 2  were  unsuccessful. 
It  may  be  that  no  solutions  exist  which  correspond  to 
charged  particles  having  a  pure  electric  dipole  moment, 
since  in  general  the  form  of  the  solution  outside  the  charge 
distribution  is  given  by  (III-7). 


1  Chapter  V,  pp.  42-46  inclusive. 

2  A  sample  result,  for  a  =  -2.177,  b  =  2.77  and  p  =  .1043 
is  plotted  (in  cross-section  for  some  fixed  values  of  r) 
in  Appendix  R. 
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VIII.  INTEGRAL  RELATIONS  SATISFIED  BY  WELL-BEHAVED  SOLUTIONS 


In  previous  chapters  it  was  shown  that  well  behaved* 
integrable  solutions  of  (V-l)  satisfy  various  integral 
relations  such  as  (V-4)*  (V-6)  and  (V-7) .  These  integral 
relations  were  obtained  with  the  aid  of  the  Lagrangian 
(V-2)  which  is  extremised  by  the  solutions  of  (V-l) .  It  is 
apparent*  however*  that  these  (and  other)  properties  of  the 
solutions  should  be  obtainable  from  (V-l)  without  recourse 
to  a  Lagrangian.  This  indeed  is  the  case. 

The  procedure  for  obtaining  integral  relations 
which  are  satisfied  by  solutions  of  (V-l)  consists*  in 
essence*  of  integrating  partially  the  identity1 * 

dn  n 

/  4>£  n  V~n  (A,<t)d35  =  /  4>|  n  -L-  F'(4>)d3|  (VII -1) 

V  1  15  V  1  n 

for  i  =  1*2*3  and  n  =  1*  2*  . . .  .  In  (VTII-l) * 

is  a  set  of  coordinates  spanning  V  and  A  ^  is  the 
Laplacian  operator  expressed  in  these  coordinates. 

1  The  use  of  this  identity  (which  follows  from  (V-l))  is. 

suggested  by  the  work  of  Morgan  and  Ladsberg  on  hypervirial 
theorems  (20) . 

The  identity  holds  for  any  closed  region  V  of  space  in 
which  the  solutions  exist  and  are  integrable  (in  the  sense 
that  all  indicated  integrals  exist).  Usually  V  will  be 
taken  to  be  all  space. 
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In  this  chapter  we  derive  some  of  the  integral 
relations  which  are  obtainable  when  the  are  taken  to 

be  rectangular  cartesian  and  spherical  polar  coordinates. 
First  of  all  we  note  that  all  integral  relations  obtainable 
from  (VTII-l)  with  n  =  0,1,2, ...,k  are  also  obtainable  from 
( VIII-l)  with  n  =  k  +  1,  since  the  "order"  n  of  the 
relation  may  be  reduced  by  one  by  partially  integrating  both 
sides  of  (VIII-l)  once.  For  the  sake  of  simplicity,  however, 
we  start  with  n  =  0:  In  this  case  (VTII-l)  implies,  simply 

/  Ab  b  d3x  =  f  b  F1  (b)  d3x  , 

V  V 

or 

/[V-(bVb)  -  (Vb)2  ]d3x  =  /  bF:,(b)  d3x  . 


Using  the  divergence  theorem1,  we  arrive  at 


§  b  Vb-ds  -  /  (Vb)2d3x  =  f  b  F1 (b)d3x  (VIII-2) 
S  V  V 


1  Throughout  this  Chapter  whenever  the  integral  theorems  of 

Gauss  and  Green  are  used  it  will  be  assumed  that  the 
solutions  of  (V-l)  and  their  indicated  derivatives  are 

well  behaved,  so  that  these  integral  theorems  are  applicable. 
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where  S  is  the  closed  surface  bounding  V.  In  particular* 
if  0  or  V4>  vanish  on  S*  (VIII-2)  becomes  (V-6),  while  for 
charged  particlelike  solutions  of  (VII-l)*  (VIII-2)  reduces 
to  (VTI-12) . 

If  is  taken  to  be  x*  one  of  the  rectangular 
cartesian  coordinates  then*  for  n  =  1*  (VIII-l)  is 


or 


f  <t>x  (A(P)d^x  =J  *x  F1  (<t>)d3x  * 


/  $x  A(|^)d3x  =  -f  x  [F (4>)-4>  F1  (^)  ]  d3x. 


V 


V 


■Sx 


Using  Green *s  theorem  on  the  left  and  integrating  partially 
on  the  right1 : 


A(x<D)d3x  +  /  [x  *  V(||)  -  ||  V(x*)]  •  da 

V  s 


(viii -3) 


=  -  /  [x  (f(^)-0  F'1  ($)}]  dydz  +  /  [ F ( ^ )  -4>  F',J ( ^ )  ]  d^x 

V  ' 


x- 


1  We  assume*  for  definiteness*  that  ^  <  x  <  x2  in  V. 
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But 


jy  ||  A(x4>)d3x  =  /v||  (xA4>  +  2  |f)d3x 


=/vtx|^F(4))  +  2(||)2  ]d3x  since  A^F1  ( $ ) , 


^  cix 


=  /  [xF^)]*2  dydz  +  /[2(||)2-  P(<t>)  ]  d3x. 


X- 


after  integration  by  parts. 
This  reduces  (VIII-3)  to 


2  Vs)2  d3x  +  A[x0y(^)  -  It  v(x<>)]  •  ds 


c3x'  3x 


(vili -4) 


“2 

:1 


=  f  [x  { 4>F"*  (4>)  -  2F( 4> )  ]  ]  dy  dz  +/V[2F(4>) -4>F*  (4>)  ]  d3x  . 

X- 


Similar  expressions  hold  for  y  and  z.  For  localised1 
solutions,  (VIII-4)  reduces  to  (V-j),  while  for  charged 
particlelike  solutions  of  (VII-l)  it  becomes  (VII-ll). 


1  Throughout  this  chapter  a  "localised  solution"  is  one  which 

vanishes  rapidly  enough  as  r  -»>  00  ,  so  that  all  indicated 
surface  integrals  vanish.  Thus  neutral  particlelike  solutions 

of  (III-4)  are  "localised"  in  this  sense,  while  the 
charged  particlelike  solutions  (p.  65  )  are  not. 
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For  n  =  2,  (VIII-l)  takes  on  the  form 

2  2 
0x2  — (A0)  d^x  =  /  $  x2  — £  F*  (0)  d^x 

8x  8x 


or 

/  <fx2A  (^|)  d3x  =  /  0>x2  [f'"  (<t’)(|^)2+  F1' 

dx 

Using  Green’s  theorem  on  the  left,  we  get: 

/v  Jj-  (x2A0+  4x4|  +  20) d^x  +  $s  [  0x2V (— p)  ^-V ( 0x£ 

Sx  8x  8x 

=  ;v  [4>X2  F1  "(*)  (||)2  +  F1 '  (■)>)  ]  d3x  • 


Since  A0  =  F1 (0)  this  becomes 


V2x  iU(^)2]  +  24,|y  ]d3x  +  ^s['x2  v 


=  /  X2[  (0F41  -  Fj  )  +  F'm.  (||-)20  ]  d3x 

8x 


or 


2 

/  x2  [0  F1  -  2F]  d3x 


8x 


( 0  )~—2'  ^  d^x 
8x 


)  ]  •  diT 


( 0x2 )  ]  •  ds* 


Integrating  partially  on  the  left  (once)  and  on  the  right 
(twice)  we  arrive  at 
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2  •/*  +  dydz 

xl 


4  f  ^^d3x  +  ^  [(t)x2V(^-|-) 

S  dx 


— «•  V(dx2)  ]  •  diT 
dx 


(VIII-5) 


x2 


J  [x2  ( cjjpl  -2F)  -  2x ( ^F1  -2F)  ]  dydz  +  2  /  ( bF'1  -2P)  d3x  , 


x- 


For  localised  solutions  of  (V-l)  and  for  the  charged  particle¬ 
like  solutions  of  (VII -1),  ( VIII-5)  is  identical  to  (VIII-4) . 

Henceforth  we  consider  only  localised,  global  solutions. 
For  n  =  3s  (VIII-l)  may  be  written  as 


/  »x3A(^i)  d3x  =  /  x3  *  [(|i)3  fiv  +  3(|i)(V|)Fin+Vl  Fi.  ]d3x 

5x3  ax  ax  Sxd 

Using  Green's  theorem  on  the  left,  and  A<t>=  F*  (<t>)  we  get 


2,530,  3$  3_ 


Va  *^3. 


6  f  x  ( — =t)  (-vr)d  X  +  6  /  X — t  ^d  x  - 
c)xd  dx  Bxd 


/  x3  [<k(||)3  fiv  +  3*  f"1  +  (^f1  1  -f1  )|-4  ]  a 


2^  ‘  1  _n3 


(3X 


3 


x  . 


Integrating  partially  on  the  left  : 
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r  r  d2(^  ,  2  d2b  ,  0  d<N  3 

■6  J  ij  (X  g-  +  2x  -^)  d  X 


dx 


dx 


,  rd2(J>  r  a,  ,  d<K  ,3 
6  J—£  (0  +  X  d  x 

dx 


-z-3 


x~  [-A  (0  F1  -  2P) 
dx3 


-  (§|)3  F1"]  d3x  , 


or 


6  /  x2  (^-|-)2  d3x  +9  / 
dx 


lx  [l  (H>2  ]d3*  +  6  ^  *3 


X 


6  /  ( ^F'1  -  2F)d3x  +  /  x3  (-^)3  F41*  d3x  • 


Further  partial  integration  of  the  integrals  on  the  left 
implies  that 


2  (^-4)2  d3x  -15  /  (l|)2  d3x  =  6/  (4>F4 -  2F)d3x  + 


6  /  x  ^ — £ 

dx 


/  X3(||)3  F>«  d3x  . 


Finally,  using  (VTII-4),  we  obtain  the  relation 


/(||)2  d3x  =  2  /  [x2  (^)2  -  \  x3  (|^)3  F1"  ]  d3x  ,  (VIII-6) 


dx 


z  A  W 


and  similar  expressions  for  y  and  z. 
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We  consider  now  integral  relations  satisfied 
by  localised,  global  solutions  of  (V-l),  when  4^  is  taken 
to  be  r,  the  radial  spherical  polar  coordinate.  When 
n  =  0,  we  obtain  the  relation  (V-6),  and  when  n  =  1  or  2, 
the  relation  (V-4)  which  also  follows  from  ( VIII -4) .  In 
any  case  no  new  relations  are  obtained  when  n  =  0,  1,  2. 


For  n  =  3^  (VIII-l)  becomes 


/  4>r3  ^  (A4>)  dJx  =  /  <t>r5  ^3  P*(*)  d°x 


(VIII-7) 


But 


,  (A*)  =  A(^i)  -  |  A  ($!*)  +  if  A(|i)  -  §  A* 
3r3  r  Sr  dr  rd 

6  S4<t>  12  33<t>  ^  36  a*  1 

+  FTTT'_?r3'  +  _zr3?  • 

or  r  or  r 


Making  use  of  Green*  s  theorem,  we  find  that 


/  br^  A (“4")  c^x  =  I  ^“3  (r^Ab  +  6r2  +  12  rb)  d^x, 

dr"5  6r 

/  br2  A(— ^-)  d^x  =  /  — g  (  r2Ab  +  4r  ^  +  6b)  d^x  , 
dr  dr 


1  This  is  shown  in  Appendix  S. 
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and 


f  brA  (||)d3x  =  /||  (rA4>  +  2  ||  +  |  0)  d~ 


x 


Thus,  ( VIII -7)  becomes 


/  C  6r 


2  S^  S<}>  oil  S2b  S$ 

5?®-24ri7® 

6  r2*  ^  + 

Sr 


36  ^  +  36  Ai)2  +  l2  + 

^2  ^  vory  r  or 


Cr3  -  6r2  ^  +  I8r  ||  -  240  J&0  ]d3x  = 


r  — V 
Sr~^ 


Sr 


/  r3  *  [(|i)3  piv  +  3  ||  ^0  pin  +  pH  ]  d3x  , 


Sr 


Sr' 


or 


r  vs  2  S3d  Sb 
j  L6r  — JBF 


24r 


S24>  Sd 


S20 


Sr' 


Sr 


2"  3r 


SdN 2  .  72  ^  Sb 


-  3 60  —2  +  36(-gjr)  +  —  T 

or 


^  r  2.  8% 
+  6  r  0  — +r 

Sr 


-  6(r 


2^  £  -3r  +  4<i>)  Ab  ]d3x  (VIII-8) 


Sr 


=  /  r3  [0  (|i)3  pIV  +  3*  ^0  p,„  +  (,p4,  _p))5V  ]  d3x  . 


^^7 


}7? 
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Integrating  by  parts : 


/-2<§S+*  f>3* 


1  3F<*§)a31 


■4  /  H>  d3x. 

It3 


3* •  °  3  /*  32o,„3..  _  ,  r  *S2(,>  „3 


/r(^|i+^)d3x=/r|_(^ 


a7  ^  T  1? 


dr 


dr' 


d  x 


and 


/  I1  d3x 
J  r  or 


\  /'*’  §p  d(r2)d« 


-  1/  [^)2  +  *f4  3 


where  r2drd^  =  d3x  .  Thus  (VIII -8)  reduces  to 


-6  /  [r2  -  3r^  +  4^]Ad  d3x  =  /  r3[ 

dr 


d3 


dF*  -2P)  -r3  (•g”) ]d3x. 


dr 


Using  Ad  =  f* (d)  and  integrating  partially  on  the  right , 


J  [r2  _  3r|i  +  4d]  F1  (d)  d3x  =  /  [10 ( dF1 -2F)+^  r3(|^)3FiM  ]d3x. 

dr 


1  J3 /ddv 3  hi  n-3. 


Finally,  using  (VIII-4),  we  obtain  the  relation: 


i  f  (v<t>)2d3x=  lo  /[  |  r3(H)3F"'  -(r2  -^-3r|^+  4f)F' ]d3x 


(VXII-9) 
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For  the  particular  case  F*  (4>)  =  $  -  (V-4),  (V-6) 

and  some  further  partial  integration  reduces  ( VIII -9)  to: 

i  /  (V4>)1 2d3x  =  |  /[<t(|£)3  r3  +  r2(3<f2-l)  (|p)2  ]d3x  .(VIII -10) 

Still  other  integral  relations  may  be  obtained  from  (VIII-l) 
by  using  higher  values  of  n  and  other  coordinates  . 

Some  uses  of  the  integral  relations  have  already 
been  demonstrated  in  previous  chapters.  They  may  be  used, 
for  example,  to  test  approximate  solutions  of  Ab  =  F'*(4>). 

o 

Thus  for  neutral  particlelike  solutions  of  A  0  =  0  -  0  , 
the  integral  relations  (V-4),  (V-6),  (VTII-6),  (V-7)  and 
(VIII -10)  imply  that  the  energy  of  the  system  represented 
by  a  solution  4> : 


M  =  ^  I  [|(V<1>)2+  -1  d]  d3x 

may  be  equivalently  expressed  by 

e  =  /<t’2  d3x  =  i b  J'<t>4  d3x  =  1 b^^)2  d3x> 

Ex  =  Tpjj:  (^r)  d^x  and  similarly  for  y  and  z. 


1  Note  that  (VIII -10)  does  not  follow  from  (VIII-6),  but 

is  (in  general)  an  independent  relation. 


' 


■ 


/ 


' 


■ 

' 


xb  ec  s  .*  a :  q>v1:>  240  tvs-iq  'r  I  .  ^^:.>:oF9b  -need 

. 

■ 

, 

. 


L 


'< 


1 

..■  ■• '  •: 


» -  - 


^ndfinaqebni  ru>  (Xjsionag  nl)  at 


96 


Er  =  5?  I  /l>  (^)3  r3  +  r2(3*2  -l)(|p)2  U3x  , 

J  or 


or 


Exx  =  h  !  [x2  +  *3  0  (H) 3  ]d3x 


and  similar  expressions  for  y  and  z. 


In  table  VIII  A,  we  list  some  of  these  quantities 
for  various  variational  approximations  to  the  two  lowest 
(in  energy)  spherically  symmetric  neutral  particlelike  solu¬ 
tions.  These  variational  approximations  were  obtained  by 
Betts,  Schiff  and  Strickfadden  (10) . 

In  table  VIII  B  we  list  similar  results  for  various 
variational  approximations  to  possible  axially  symmetric, 
odd  parity  neutral  particlelike  eigensolutions .  The  indicated 
approximate  variational  solutions  are  those  which  were 
listed  previously  in  tables  VI  A  and  VI  B. 

mc^tt 

In  these  tables  all  results  are  in  units  of  — . 

6 

The  results  for  the  nonspherical  functions  were  obtained 
in  part  by  numerical  integration.  The  symbol  E,  where  it 
appears,  implies  that  the  corresponding  integral  relation  is 
identically  satisfied  by  the  indicated  trial  function.  For 
example,  in  table  VTII  B,  E  appears  under  E^  and  E^ 
for  the  function  Az T  e~rT.  This  means,  simply,  that  this 
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trial  function  satisfies  identically  the  integral  relation 

E  =  E  =  E 
x  z 

It  is  interesting  to  note,  in  table  VIII  A,  that 

if  AE  represents  the  difference  (in  absolute  value) 

between  the  values  of  E  of  a  given  trial  function  and  the 

"exact”  solution,  then  the  corresponding  deviations  for  the 

other  integrals  (that  is  AE  ,  AE  ,  etc.)  are  at  most  one 

x  z 

order  greater  than  AE.  For  example  for  the  trial  function  2 
AE  ~  .04  while  AE  ~  .12  and  AE  ~  .08  and  for  the 

~  ±  XX 

function  5*  AE  ~  .51  while  /AE  ~  1.36  and  AE  =  1.46. 

—  X*  XX 

This  is  consistent  with  the  fact  that  the  trial  functions 
extremise  E  but  not  (necessarily)  E  ,  E  ,  etc. 

X  Zj 

It  is  not  clear  whether  the  results  for  the  odd  parity 

'  -  f-~)  Y*  -•  '  Y  f  ’!  )  *  *  •  *  •  -  '  :  ; 

states  exhibit  such  a  consistency.  As  pointed  out  in  Chapter  VI 

some  of  these  variational  trial  functions,  especially  4>B, 

and  seem  to  approximate  a  two-particle  state  rather  than  a 

ill 

one-particle  state  of  odd  parity.  It  is  clear,  also,  that  non- 
spherical  solutions  are  much  more  difficult  to  approximate 
vari ati onally  than  spherical  ones,  hence  these  approximate 
results  may  still  be  appreciably  different  from  the  exact 
solutions.  As  a  result,  a  trial  function,  which  gives  a  reasonable 
approximation  to  the  energy  E,  may  still  give  poor  results  for 
the  other  integrals,  which  are  not  extremised  in  the  variation. 
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IX.  STABILITY  OF  PARTICLE LIKE  SOLUTIONS 


The  question  of  stability  of  particlelike  solutions 


of 


Ab  =  F  (b) 


(IX-1) 


has  been  considered  by  Derrick  (l8)  who  proved  (as  already 
mentioned  in  Chapter  V)  that  the  integrable  solutions  are 

unstable  with  respect  to  radial  perturbations.  Rosen  (21 ) 

!  1 

then  proved  that  Derrick  s  instability  criterion  implies,, 

2 

for  a  wide  class  of  field  equations.,  that  the  solutions  are 
also  dynamically  unstable  in  the  more  general  sense  of 
Liapunov. 


In  this  chapter  we  consider  the  dynamical  stability 
of  particlelike  solutions  of  the  field  equations  (i-l).  More 
specifically  we  shall  consider  only  the  dynamical  stability 
of  the  neutral  particlelike  solutions  discussed  in  Chapter  VI. 
These  solutions  are  not  stable  in  the  sense  of  Derrick  s  theorem. 


1  vi s  8  L  [ b  (pr) ] 


p=l 


<  0  for  all  integrable  solutions 

b (r)  . 


2  i.e.,  those  derivable  from  the  Lagrangian  density 


X  =  (v<t>)2  -  (|£)2  + (*). 
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however  their  (physically  more  meaningful)  dynamical  stability 
is  not  immediately  precluded.  This  is  because  the  Lagrangian 
(1-2)  does  not  fall  within  the  class  =  (V4>)  ■(^t) +F  (<l>) 
mentioned  in  the  footnote  on  p.  101 3  even  when  if  =  0.  ^ 

Actually,  one  might  expect  the  neutral  particle¬ 
like  solutions  to  be  unstable,,  particularly  so  since,  as 
pointed  out  by  Rosen  (21),  all  spinless  "elementary"  particles 

(such  as  t r,  K)  known  to  date  are  only  "metastable”  (that  is, 

2 

their  decay  rate  is  «  . 

To  determine  the  dynamical  stability  of  a  neutral 
particlelike  solution  <f>Q(r),  we  consider  a  general  time- 
dependent  perturbation  of  this  solution: 


A±  =  [cf,  i  <t>0 (r)  ]  ->  A±  +  5A±  =  [a(r,t),  i  {dQ(r)+^(r,t) }  ]  . 


We  assume  that  the  perturbations  a  and  ^  are  initially  small. 
Substituting  the  perturbed  expressions  into  the  field  equations 
(i-la),  we  get  (to  first  order  of  small  quantities): 


1 


In  fact 
A  =  0, 


,  the  field  equations  (i-l) 
0,  that  is  only  static 


imply  that  when 
solutions  are  possible. 
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V  x  V  x  a  + 


+ 


+  1^-  (V-a)  =  (1  -  3*02)^ 


(IX-2) 


From  physical  considerations  we  expect  the  particles  to 
decay  into  other  particles  or/and  radiation. 

Consider  first  the  possibility  that  the  particle 
decays  into  (plane  wave)  radiation  directly.  Thus.,  we  assume 
that  the  (initially  small)  perturbation  is  plane  wave-like 
only.1  In  general.,  then: 


^  =  S  o  ak  (ak) 
k 

and 

*  =  f  ck  he 

k 

where  ak  =  £*r  -  ^t  and  the  ck  are  arbitrary  but  small 
constants.  The  summations  are  taken  over  all  values  of  the 
four-vector  k  =  (It,  i^)  . 

2 

Substituting  these  expressions  into  (IX-2)  we  get  , 
by  virtue  of  the  independence  of  the  constants  ck, 

1  Note  that  this  is  not  a  localised  perturbation. 

2  Recall  that  V  x  V  x  a  =  V(V’a)  -  Aa  . 
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f5*V)V  sc  7  x  7  x  7  lartf  IlBoefl'  S 


io4 

ic(Ic«ak  )  -  k2ak  -  ^k  -  (<t>02-l)ak  =  0  (IX"3) 

and 

ic2^  ~  CD^’*ak  +  (3  <t>o2_1^k  =  (IX-4) 

d  2  -*2.  2 

where  the  prime  —  and  k  =  k  -  a)  .  Since  the 

k 

coefficient  of  ^k  in  (IX-4)  cannot  he  expressed  as  a 
function  of  ak  only,  we  must  have  ^k  =  0,  hence  (for 

Tf 

nontrivial  solutions)  lc-ak  =  0.  (IX-3)  then  reduce  to 

k2  ak"  +  [<t>02(r>)  -  1]  ak  =  0 
which  again  contradicts  our  assumption  that 

a,  =  a,  (cr.  )  only.  Thus,  the  initial  perturbation  cannot  be 
K!  xC 

purely  plane  wave -like,  implying  that  it  is  impossible  for 
a  neutral  particle  to  decay  directly  into  plane  wave  radiation. 

Consider  now  the  possibility  that  the  particle 

represented  by  ^(r)  decays  directly  into  another  particle. 

Assume  that  the  perturbations  can  be  expanded  in  a  series 

i  kt 

over  the  (complete)  set  of  functions  e  : 

a  =  2  Re  uk(r)  ckeikt  , 

and  (IX-5) 

ip  =  2  Re  v1_(r)ck eikt 

k  k  K 


JbnB 


..  ©fffinq  9dcf  WlSiivt 


?a*jm  ew  '%X»no-  .  p  'io 


noivtonjj'i 


.  0  ~  *6*1  'ir'oxn ij  i  ;  i&tvi-.  3 non 


, 


■ 


,L  : 


: 


o.  :zc  1 1  ‘  .  .lv  j;  .  .  jXnd1"'  -L{o}  t/  •-  £ 


3  ;  ;r.  y  :*-.r i \<£qe! t  tS3j£I-$fVj$W  sfTSlq 


,  i  can,  ,1  'o fVT  Lb  ’OS  9  Xtajsq  XirrtJJdO  £ 


aXS-tJ' ijsq  srij  ^ifilXcfieeoq  erlcr  won  ^dfclenoO 

V 

>  £■  is: 0  3b  ( *i )  ( .. ♦  <J5  ;  j n9s erxqsn 


enol.'onjjl  lo  368  vetfalqmoa)  ari3  lavo 


\ 


105 


where  the  c^  are  arbitrary  (small)  constants.  The  summa¬ 
tion  is  taken  over  all  admissible  values  of  k.  Substituting 
(IX-5)  into  ( IX-2) ,  we  find  that  uk  and  vR  must  satisfy: 


V  x  V  x  uk  -  (k2  +  4>02  -  l)uk  +  ik  V  vk  =  0 

o  (ix-6) 

A  vk  +  *0  '  1^Vk  +  lk  V‘"k  =  0  • 


Since  the  perturbations  are  particlelike,  we  require  that 

u^  and  v^  be  everywhere  finite,1  and  u^,  v^  must  ->  0 

2 

as  r  _  oo  for  a  localised  disturbance. 

The  problem  is  thus  to  find  the  eigenvalues  k  of 
this  set  of  coupled  equations  (IX-6),  with  the  indicated 
boundary  conditions.  An  imaginary  eigenvalue,  if  one  turned 
up,  would  imply  a  decay  rate  of  (-k2)1/2  . 

In  general  it  appears  rather  difficult  to  say 
anything  definite  about  the  nature  of  the  eigenvalue  spectrum 
of  k  (e.g.  whether  discrete,  real,  etc.).  We  consider  some 


1  Also  since  the  perturbations  are  assumed  initially  small. 

2  This  might  be  relaxed  somewhat  to  requiring  that  vk  -*  a  constant 
yk  (for  example  for  a  charged  particlelike  perturbation). 
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special  forms  of  the  disturbance  which  reduce  (IX-6)  to  a 
simpler  form.  Assume  a  perturbation  for  which  V  x  u^  =  0 
that  is  H  remains  null.  It  then  follows  from  (IX-6)  that 


ik  Vv 


uk  = 


k 


(k2-!-^2) 


and 


Avk  [1 


k^ 


2  2 
k  -i+*0 


] 


2k24>nO  3vk  p 

+  s/  +(3^02-l)vk 


=  0 


(k  -l+4>0  ) 


(IX-7) 


Let  vk  =  S  Ifon  wk(1')Yi;(0>'!))  >  being  arbitrary 

&  ,m 

constants.  Then,  from  (lX-7).>  we  obtain 


1. 


d2w. 


dr 


|+2  [i  + 


k2^  0* 

0  0 

(4>02-l)(k2+<t>02-l) 


] 


dw. 


k 


dr 


(3b0  )(k  -1+*Q  )  t(l+l) 

+  L - 2 - I - 2^Jwk 

0O  -  1 


=  0  (ix-7a) 

with  w^(r)  0  as  r  -*  00 

and  lim  C  r  w^(r)]  =  0  . 

r-+0 


Approximate  solutions  of  this  equation  were  obtained 

2 

numerically.  One  such  eigensolution  (with  k  =  -14.53  and 

i  =  0)  is  plotted  in  figure  9.1.  No  well  behaved  eigensolutions 
2 

were  found,  thus  contradicting  our  initial  assumption  of 


1 

2 


,  d0 
0’  =  __ 2. 
0  dr. 


Note  that  the  coefficients  of  the  equation  are  singular  at 
o 

0  =  1.  The  eigensolutions  also  appear  to  be  singular  there. 
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possible  it  ?=  0,  initially  small  perturbations.  It  is  not 
possible,  therefore,  for  neutral  particles  to  decay  directly 
into  particles  for  which  it  0,  that  is  they  are  stable 
against  this  decay  mode. 

Further  inspection  of  (IX-6)  indicates  also  that 
neutral  particles  are  stable  against  a  decay  mode  in  which 
either  the  vector  or  scalar  potentials  remain  unperturbed 
(i  .e  .  uk  or  vk  =:  0)  . 

We  consider,  finally,  a  perturbation  for  which 
V*uk  =  0,  whereupon  (IX-6)  reduce  to 

Affk  +  (k2  -  1  +  002)uk  +  ik  V  vk  =  0 

and 

Avk  +  (3*02  -  1)  vk  -  0 

The  latter  equation,  with  the  boundary  conditions 

lim  v,  (r)  =0  and  lim  [r  vk(r) ]  =  0  has  only  the  trivial 
r->  “  K  r^O 

solution  vk  =  0,  in  which  case  (IX-8)  becomes 

Auk  +  (k2  -  1  +  >f02)uk  =  0  (IX-10) 

with  lim  u,  (r)  =  0  ,  lim  r  u.  (r)  =  0 

r-*  oo  K  r-»0 

and  V*uk  = 


(IX-8) 

(IX-9) 


0 


(IX-11) 
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A  solution  of  (iX-ll)  is  uk(r)  =  uk(r,0)  ^  -> 

Tj  being  a  unit  vector  in  the  direction  of  increasing  cp. 

cp 

(lX-10)  then  becomes 

Auk  +  (K2  +  0Q2  -  ~  *  2  )uk  =  0  (IX-12) 

r  sm  0 

2  2 

where  K  =  k  -  1  .  ( IX-12)  can  be  simplified  by  variable 

separation,  whereupon  we  find  that 


00 


uk(r,0)  =  2  -  fk<e(r)  (cos  0)  , 

£— 1 


£l 


where  P^  are  the  associated  Legendre  functions  P 
with  m  =  1,  and  fk^(r)  are  eigensolutions  of 


£m 


d2f 


k  £ 


dr 


+  [k2  +  <t>Q2(r) 


l ( ^+l)  1  f 

2  J  1ki 
r 


=  0 


f 


with  i  >  1,  f (0)  =  0  and  lim  (-—■)  =  0  . 


p-,^0  r 

When  i  =  0,  the  equation  (IX-13)  has  the  solution 

.2 


(IX-13) 


f 


.  =  r  <t>  and  YT  =  -1  hence  the  ground  state  eigensolution 

ko  0 

2  2 

of  (IX-13)  (when  £  >  l)  is  clearly  >-  1,  that  is  k  >  0. 

It  then  follows  that  the  (spherically  symmetric)  neutral 
particles  are  stable  against  this  particular  decay  mode. 


As  remarked  previously,  it  is  not  essential  to 
consider  localised  particlelike  perturbations  only.  Thus  the 
boundary  condition  on  (IX-6)  a t  infinity  can  be  generalised  to 
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uk,  vk  -►  const.  /  0  as  r  —  oo  ^ 

corresponding  to  a  nonlocalised,  but  still  particlelike 
and  initially  small  (by  virtue  of  the  small  coefficients 
ck  which  appear  in  the  expansions  (IX-5)  )  perturbation. 

This,  however,  reduces  to  the  localised  situation,  at  least 
for  the  two  special  cases  of  ()IX-6)  which  have  been  considered 
(V  x  u,  =0  and  V«u.  =  0) .  This  is  because  for  large  r 
(when  $0  0  as  e~r/r )  (IX-7)  and  (IX-9)  reduce  to 


A  vk  +  (k2  -  1)  vk  =  0 

and 


A  vk  -  vk  =  0  respectively,  and 
the  well  behaved  solutions  of  these  equations  are  asymptotic 


to  zero  at  infinity. 
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X.  CONCLUSIONS 

In  this  thesis  we  have  been  concerned,  primarily., 
with  the  investigation  of  the  properties  of  static,,  particlelike 
solutions  (in  the  absence  of  the  magnetic  field)  of  the  field 
equations  proposed  by  Dr.  Schiff.  The  only  other  types  of 
solutions  investigated  at  some  length  were  plane -wavelike 
solutions.  It  was  found  that,  in  addition  to  the  usual  Maxwellian 
plane  waves,  the  theory  admits  plane -wave like  solutions  for 
which  the  electric  field  is  not  orthogonal  to  the  direction  of 
propagation  of  the  wave  and  the  charge  and  current  densities 
are  non-zero. 

The  existence  of  global,  well  behaved,  spherically 

Q 

symmetrical  solutions  of  =  <$>  -  was  formally  established. 

The  existence  of  well  behaved  nonspherical  solutions,  however, 
was  demonstrated  only  in  the  vicinity  of  an  arbitrary  point  of 
space . 

The  variational  method  for  obtaining  approximations 
to  well  behaved  solutions  of  the  general  class  of  field 
equations  A0  =  FT(b)  was  considered  and  sufficient  conditions 
were  derived  under  which  a  (variational)  upper  bound  to  the 
Lagrangian  is  assured.  These  conditions  assure  upper  bounds 
for  all  eigenstates  (  not  only  the  one  of  lowest  energy  ) 
without  recoure  to  any  other  subsidiary  conditions  (  such  as 
orthogonality  conditions).  The  usual  technique  of  parameter 
variation  was  generalised,  for  the  non-spherical  case,  to 
include  variation  of  function  of  one  of  the  independent  variables. 
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The  list  of  known  spherically  symmetric  neutral  and 
charged  particlelike  solutions  of  the  field  equations  was 
extended  to  include  ’'compound"  particlelike  solutions,  which 
are  linear  in  the  vicinity  of  (some  of)  their  zeros. 

A  variational  principle  for  the  charged  particlelike 
solutions  was  determined  and  used  to  obtain  approximations  to 
the  lowest  (in  energy)  spherically  symmetric  charged  particlelike 
state.  It  was  found  that  this  procedure  is  rather  inefficient 
for  obtaining  approximations  to  the  mass  and  charge  of  the 
particle,  since  it  is  the  Lagrangian  and  not  these  quantities 
which  is  extremised  in  the  variational  procedure. 


Approximation  techniques,  numerical  and  variational, 

were  used  in  an  attempt  to  determine  non-spherical,  axially 

symmetric,  neutral  particlelike  solutions  of  odd  parity.  The 

results,  unfortunately,  do  not  unambiguously  indicate  the 

existence  of  neutral,  one-particle,  non-spherical  states.  The 

numerical  results  (obtained  by  iteratively  solving  the  finite 

difference  approximation  to  the  field  equation)  could  not,  in 

any  sense,  be  interpreted  as  representing  well  behaved  solutions, 

since  they  were  found  to  be  highly  unstable  near  the  origin. 

The  variational  results  on  the  other  hand,  seem  to  suggest  a 

number  of  possibilities:  The  best  results  obtained  by  the 

parameter  variation  technique  (table  VI -A)  suggest  the  existence 

me 

of  a  one -particle,  odd  parity  state  of  mass  <  3-9 6  — ? —  . 
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The  nodeless  variational  trial  functions,  obtained  by  the 

variation  of  function  method,  (table  VI-B),  it  is  believed, 

approximate  the  lowest  two-particle  state  which  has  an  energy 
,  me 

of  3*006  — — .  It  is  possible,  however,  that  they  approximate 

g 

a  one  particle  odd  parity  state  of  energy  not  much  above  this 
value.  The  various  trial  functions  with  nodes  (table  VT-B), 
are  approximations,  in  all  probability,  to  various  excited 
odd  parity  states. 


A  numerical  search  was  also  made  for  a  (non-spherical ) 
particlelike  state  corresponding  to  a  charged  particle  with 
a  pure  electric  dipole  moment.  No  such  solution  was  found. 


A  scheme  for  obtaining  integral  relations  satisfied 
by  integrable  solutions  of  AO  —  FT(0)  was  developed  and  used 
to  obtain  the  first  few  integral  relations.  These  relations 
were  then  used  to  obtain  various  expressions  for  the  energy 
corresponding  to  a  given  solution.  They  were  used,  also,  to 

3 

test  variational  approximations  to  the  solutions  of  AO  =  0-0  . 
As  might  be  expected,  it  was  found  that  the  variational 
approximations  to  the  spherical  states  satisfy  these  integral 
relations  much  more  closely  than  do  approximations  to  odd 
parity  states. 

Finally,  the  dynamical  stability  of  the  spherically 
symmetric  particlelike  solutions  was  investigated.  They  were 
found  to  be  stable  against  many  dissociation  modes,  among  them 
against  dissociation  into  plane  waves  and  particles  with  H  =  0. 
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APPENDIX  A:  COMPONENTS  OP  ENERGY -MOMENTUM  TENSOR 


In  vectorial  notation,  the  components  of  the 
symmetrised  energy -momentum  tensor  (1-3)*  are: 


T„3  -  fa  [?(e2+s2)6aP-  EaEp  '  HaH@  ‘  +  2 ) 2+K )S 

+  (1  +  S2  -  P)AaAp  ] 

TUa  =  fa  C  (ff  x  ff)a  +  (1  +  S2  -  *2)A0*]  , 

t44  =  -  ^  +  if )  +  |(^£-‘t'2)  +  ^(if-o2)2  +  K 

+  (1  +  S2  -  4)2)4>2]  . 
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APPENDIX  B:  RELATION  BETWEEN  REST-MASS  AND  LAGRANGIAN 


The  energy  of  a  system  characterised  by  the 
potentials  is  (in  the  rest  frame): 
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M  =  -  /T44d  x  =  7*7 


r  r  1  d  2  1A  2  .  1  A  2A  2 

J  Llf  +  “A-*  +  ^  Aj  A’ 
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¥  Hi 


(1  +  Ak2)A42  +  K]  d3x 


From  the  field  equations,  (i-l) ,  we  obtain 

Q  Q  p  "5 

f  A1Fik  k  d^x  =-/(!  +  Ak‘~)Ai  dJx  (no  summation  on  i )  . 
Since  AiPlkjk  =  (Vile*,*  -  Ai,kPik  and 
/  (Aipikbkd3x  =  /  (AiPiaba  d3x  +  *  (AiPi4b4  d3x 

=  0  Ai Pi ads a  +  /  (AiPiV,4  ^  ' 


the  above  becomes 


0  AiPiadsa  +  A.(Vl4>,4-  Ai,kPik  +  +  Ak2)A.2]d3x  =  0  . 
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Equivalently : 
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Ak,iFik 
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With  i  =  4,  this  reduces  the  expression  for  the  rest  mass  to 

Fik2  +  I  Ai2  +  ¥  AiV  +  K^x  ' 

'  ¥7  /  Aa,4F4ad3x  +  ¥?  ^  A4F4adsa  ' 
which  is  the  relation  (1-12). 


1.  For  the  static  case,  these  integral  relations  follow 
from  (1-10) . 
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APPENDIX  C:  SOME  UNPHYSICAL  SOLUTIONS 


The  discussion  of  plane -wave  solutions  indicates 

the  existence  of  analogous  "standing  wave"  and  "homogeneous" 

— ► 

solutions  corresponding  to  ^  =  0  and  k  =  0  respectively. 

Many  of  these  solutions  are  just  limiting  cases  of  the 
plane  wave  solutions.  A  large  number  of  unphysical  solutions 
may  also  be  obtained  by  using  the  usual  cylindrical  coordinates 
(p,  <p,  z).  For  example,  if  1  =  0,  the  field  equations 
reduce  to: 


id  /  d<l>\  \ _  d^4>  ,  d^4> 

p  ^  (P^}  +  +  ^5 


-  <J)  +  =  0 


P  dcp* 


dz‘ 


/  d  ^  d  ^  \ 

For  the  case  of  cylindrical  symmetry  (when  ^  ^  =  0) 

at  least,  phase  space  analysis  of  the  type  used  by  Finkelstein 
et  al  (15) »  indicates  the  existence  of  a  countable  (discrete) 


set  of  well  behaved  solutions  asymptotic  to  zero  as 
and  a  continuum  of  solutions  asymptotic  to  -  1  as 
sin(^2  p  +  c) 


e  ~P 
v'p 


Another  example  corresponds  to  =  AeT  and 
0  =  yA  where  ?  is  a  unit  vector  fixed  in  space  and  y 
a  constant.  For  the  steady  state  case  =  0),  the  field 
equations  (i-l)  become 
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eAA  -  V(e*VA)  =  (l  -  A2A2)Ae 
AA  =  (1  -  A2A2)A 


2  2 

with  A  =7  -  1. 


Choosing  the  z-axis  along  e  ,  this  implies  that 


$>A 


=  fz  =  0  and  A  =  ^(P^)  where  f  is  a  solution  of 


I  (p  |£)  +  1  aff  =  f  .  A2f3  , 


which  has  well  behaved  global  solutions  (at  least  for 
A2  >  0) . 
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APPENDIX  D:  FIELD  EQUATIONS  IN  SPHERICAL  POLAR  COORDINATES 


The  field  equations  (i-l),  in  spherical  polar 
coordinates. 


x  =  :(r,0,9) 


^  (Ar,A0,Acp) 


take  on  the  form: 
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APPENDIX  E;  UNIQUENESS  OF  THE  SPHERICALLY  SYMMETRIC  CAUCHY 

PROBLEM  FOR  A0  =  0  -  03 


Consider  the  Cauchy  problem: 


Ab  =  0  -  for  r  <  R 


with  0 (R)  =  -  1, 


^0 

^r 


r=R 


=  dj  R  >  0  and  d  being  constants. 


Let  0(r)  =  -  1  +  x(r)  be  a  solution  of  the  problem.  X(r) 
then  satisfies 


AX  +  2x  -  X2  +  X3  =  0 

with  x  =  °,  ^  =  d  on  r  =  R. 

Near  the  surface  r  =  R  (where  |  x  I  «  1  ) 

00  £ 

X(r)  -  S  2  [Aim  J^(x)  +  B^n  (x)  ]  y“  (6,9), 

£=0  m=-£ 

where  x  =  ,/2  v3  and  n^  are  spherical  Bessel  and  Neumann 
functions  of  order  £  respectively  (Ref.  (17),  p.  539),  and 
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are  constants. 
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At  r  =  R,  x(R->0j9) 


?n  j  .  [klmhM  +  BJmn  (x)]Y“(0j9)  =  0 

£=0  m=-i> 


hence  A^m^(x)  +  Bxmni(x)  =  0  for  a11  ^  and  m> 

where  X  =  y/2  R  . 


Similarly, 


from 


d<t> 

3r 


r=R 


d,  we  get 


AtoJi  (x)  +  Bimnl  (X)  =  0  for  a11  £  >  0 

and  |  m  I  <  i.  (The  prime  denotes  differentiation  with 

respect  to  the  argument).  If  X  is  a  zero  of  one  of 

n^  or  their  first  derivatives,  then  these  conditions  imply 

that 


A,  =  B . 
im  Hm 


=  0  for  &  >  0,  |  m  |  <  & 


Otherwise 
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'Jim 


J -(X) 

fvi  j  hence 

n^(X)  Jim 


Eyxyw  [nBx)>t(x)]x=x  = 


W  being  the  Wronskian  of  n^  and  j^,  that  is 


W  [n^j^] 


x=X  2R 


(Ref.  (17),  P.  541). 
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It  then  follows  that  A^m  =  B^m  -  0  for  all  &  >  0,  |  m  |  <  i3 
so  that  4>(r)  Is  spherically  symmetric  in  the  neighbourhood 
of  r  =  R.  However  f(r)  =  r  4>(r)  satisfies  (in  the  limit 
as  Sr,  6|i-^0)  the  difference  equation  (V-13)*  from  which 
it  is  obvious  that  since  the  solution  is  spherically  symmetric 
at  r  =  R  and  r  =  R  -  Sr  (for  arbitrarily  small  6r  >  0), 
then  it  is  such  for  r  =  R  -  26r,  hence  by  induction  for 


all  r  <  R. 
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APPENDIX  F:  EXISTENCE  OF  SOLUTIONS  FOR  (IV-4) 


For  the  system  of  linear,  inhomogeneous  differen¬ 
tial  equations  (IV-4)  assume  that  the  solutions  for 
i  =  1,  2,  . ..,  i  are  polynomials  of  degree  <  i  for  each  i. 
Consider  now  the  equation  for  i  =  H  +  2: 

1. 

^ai+2^^  +  (^+3)  (^+2)u^+2(q)  =  (F-l) 


where 


Pi+2 


l  m 
Z  Z 
m=0  n=0 


a 

m~n 


a 

n 


1 

2  b(1P,W 

j=Q  " J  J 


(say) 


(F-2) 


00 


Assume  Q^+2(p)  =  2  a^.P^  (p)  =  Substituting 

into  ( F-l )  we  find  that. 
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^  t-j(j+l)a^j  +  (i+3)  (^+2)a^j-b^j]  Pj(d) 
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from  which  it  follows  that  (F-l)  has  the  solution 
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Furthermore *  this  solution  is  a  polynomial  of  degree  £, 
hence  is  analytic  everywhere  in  |  \i  |  <  1. 

For  £  =  0,  (IV-4)  has  the  solution 


a2(n)  =  C2P2(n)  +  ^(CQ-Co3) 

1  S 

which  is  of  the  form  (F-2)  as  assumed  (with  b2Q  =  -£-(Cq  “  ) j 

b2-^  =  0  and  b22  =  C2).  Thus.,  the  existence  of  well 
behaved  solutions  for  all  follows  by  induction. 
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APPENDIX  G:  PARITY  OF  SOLUTIONS  OF  (IV-4) 


As  shown  in  Chapter  IV  and  Appendix  F,  the 
solutions  of  (IV-4)  may  be  written  as 
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Substituting  into  ( G-3 )  we  get  a^(-p.)  =  (-l)kak(M-). 
(G-2)  holds  for  i-  0,1,2,  then  by  induction  it  holds 
all  i. 
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APPENDIX  H.  RADIUS  OF  CONVERGENCE 


Consider  the  radius  of  convergence  of  the  series 
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Since  (H-3)  is  true  for  i  =  1  (a^  =  £-(a0-aQ3)  )  *  then  by 

induction  (H-3)  holds  for  all  Z, 

Thus  we  have 
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But,  by  the  ratio  test,  the  latter  series  converges  uni¬ 
formly  to  an  analytic  function  for 
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hence  the  same  applies  to  (H-l). 
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APPENDIX  J.  SERIES  SOLUTION  OF  (lV-6) 


In  the  table  below  we  list  some  examples  of  the 
solutions  of  (IV-6)  as  evaluated  from  the  series  (IV-5) 
for  various  0(0).  K  is  the  number  of  terms  of  the  series 
required  for  the  tabulated  results  to  be  exact  to  (at  least) 
four  decimal  places.  The  last  line,  for  each  0(0) 


corresponds  to  that  value  of  r 
appears  to  diverge  (then  ^300  = 


for  which  the  series  first 

3Z°  a.  r2*+1). 
i=0  “ 


0(0) 


r 

0 

0.200 

0.50250 

0.400 

0.51002 

0.600 

0.52259 

0.800 

0.54028 

1.000 

0.56314 

1.200 

0.59119 

1.400 

0.62440 

1.600 

0 . 66262 

1.800 

0.70551 

2.000 

0.75256 

2.200 

0.80287 

2.400 

0.85521 

2.600 

0.90781 

2.800 

0.95877 

3.000 

1.00565 

3.200 

1.04575 

3.400 

1.07846 

3.600 

1.10047 

3.800 

1.11187 

4.000 

1.11571 

4.200 

1.14296 

4.400 

1.09H5 

4.6 

03OO=1 ,20^ 

K 

2 

2 

3 

3 

4 
4 

4 

5 

5 

6 
7 
7 
9 

10 

10 

10 

15 

20 

23 

23 

23 

117 

106  a  301=-6.313  105 


. 

i 

. 


oesoc.o 

OOh.Q 

aiede  o 

. 

•Cr  ..I 

OCJV.  £ 

-.cecv.0 

8.0 

G08.S 

' 


■ 


. 


(0)  =  1.25 


r 

0 

K 

0.200 

1.24535 

2 

0.400 

1.23179 

3 

o.6oo 

1.21047 

4 

0.800 

1.18311 

4 

1.000 

1.15168 

5 

1.200 

1.11833 

6 

1.400 

1.08495 

8 

1.600 

1.05323 

10 

1.800 

1.02445 

13 

2.000 

0.99961 

17 

2.200 

0 . 97937 

26 

2.400 

0.96393 

48 

2.600 

0.95326 

275 

2.8 

«3ocr4-485  1q1' 

a30_=8.484  10l4 

(0)  =  2 


r 

4> 

K 

0.200 

1 . 96086 

3 

0.400 

1 . 85300 

4 

0.600 

1 . 70002 

6 

0.800 

1 . 52888 

10 

1.000 

1.36137 

16 

1.200 

1.21093 

35 

1.4 

*30Cf3- 579  lc2 

a3C_ =7.067  1C“ 

(0)  =  5 


r 

4> 

K 

0.100 

4.80714 

4 

0.200 

4.30317 

6 

0.300 

3-64990 

10 

0.400 

■2.99075 

21 

0.500 

2.40473 

136 

30 


n 


OX  484.8 


01 

88 


2 


4 

01 


01  Pl£ . 1 ~p  jo 


3t 

;  :  v. . 

4 

■n  C-'.  i 

TWI2.X 

4 

U:-8i.  I 

0 

Os  -  .  i 

!  . ; 

■:  !'?.■•.  !. 

.•  ■_  t-'O.j: 

c?t>i!'20..£ 

YI 

,;P .  0  fa  .  . 

rcevt  -  o 

ut> .  ■) 

c!YS 

V.VcP.O 

8:  ^  -  u 

:Q£ 

- 

$ 

003.  o 


008 . 0 
000  .  X 

oc  :•’ .  i 

" 

008.  . 


00  ::.  U 
.004,2  ■ 
OOh,  S 


8.0 


080 oP.O 
ooeeB . i 
, 

r£lci£..X 

8^12.1 


002,  0 
004.0 

ooo. : 

OOS.l 


k0X  Ok  -e-oo.c^  4 .  !' 


.  ■• 


4> 


4IY08 . 4 

YX80r . 4 
0 

:TOQO.O 

4Y404.S 


001.0 
004.0 

00^.0  '  - 


■' 


J-3 


0(0)  =  10 


r 

$ 

K 

0.100 

8.56449 

7 

0.200 

5*  8666l 

23 

.3 

- 

a224=-!.350  10 

K-l 


APPENDIX  K.  VARIATIONAL  APPROXIMATION  TO  POSSIBLE  ODD -PARITY 

NEUTRAL  PARTICLE  EIGENSOLUTIONS  USING  THE  COMPARISON  FUNCTION 

0A  =  R(r)  P1(p.) 


Substitution  of  the  comparison  function  <J>A  =  R(r)P-^(|j.) 
into  the  Lagrangian  (VT-2),  yields 
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(K-l) 


The  corresponding  Euler-Lagrange  equation  is 
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In  analogy  to  (IV-6),  this  equation  has  a  discrete  set  of 
"orthogonal"  1 2  eigensolutions  asymptotic  to  zero  as 

a  2 ( i  +  =  A  K  frO 

i  r  '  rJ  i  Jj  which,  presumably,  yield  upper 


bounds  to  the  masses  of  the  particles  represented  by  odd- 
parity  eigensolutions  of  (III-4) . 


1  Orthogonal  in  the  sense  of  (VT-l) . 

2  The  first  three  values  of  the  asymptotic  amplitude  A.  are 
listed  in  table  (VI-B),  while  the  first  three  eigensolutions 
are  plotted  in  figure  6.6. 
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From  (K-2)  it  follows  that  the  solutions  satisfy 
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which  reduces  (K-l)  to 
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dr  . 


Note  that  the  eigensolutions  of  (K-2)  immediately  satisfy 
the  integral  relations  (V-4)  and  (V-6),  since  the  latter  are 
satisfied  by  all  comparison  functions  for  which  L  has  been 
extremised  with  respect  to  amplitude  and  radial  scale  parameters.1 2 


written 


Near  the  origin,  the  solutions  of  (K-2)  may  be 
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where  a^  is  arbitrary,  while 
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m=0  n=0 

for  Z  =  1,2,  ... 

1  Thus  (K-3)  is  obtainable  from  (V-6) . 

2  This  series  expansion  of  R  in  the  vicinity  of  r  =  0  is 
obtainable  in  the  same  fashion  as  (IV-5) • 
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The  first  three  eigensolutions  of  (K-2),  which  are  plotted  in 
figure  6.6,  were  obtained  numerically  by  solving  the  finite 
difference  approximation 

ijr+S)  +  f(r-6)'-  2f(r)  =  '  (l  +  2U(r)  _  3 

6^  r  D  r 

to  (K-2)  (f  =  rR  and  5  is  the  mesh  size),  using  the  method 
described  in  Chapter  V. 


The  integral  expressions  which  are  given  in  table 
VIII  B  were  obtained  numerically  from 
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APPENDIX  L.  VARIATIONAL  APPROXIMATIONS  TO  ODD-PARITY  NEUTRAL 
PARTICLE  El  GENS  OLUTI ONS  USING  THE  TRIAL  FUNCTION  4>B=  rne“arY(p) 


The  function  <l>  =  rn  e  arY(q)  reduces  (VT-2)  to' 
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Setting  v  =  1  -  \±,  (L-2)  becomes 

37[(2v  -  v2)§f]  -  c2y  +  Cl(Y3  =  0  (L-3) 


with  Y(v=l)=0 
Near  v  =  0  Y  = 
is  arbitrary,  while 


for  an  odd  parity  solution. 
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The  finite  difference  approximation  to  (L-3)  is 
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-C2Y(v)+C4Y3(v) 


=  0 

Numerical  solution  of  this  equation,  using  the  iteration 
technique  described  in  Chapter  V  indicates  the  existence  of  a 
discrete  set  of  orthogonal  (in  the  sense  of  VT-l)  odd-parity 
eigensolutions  of  (L-2)  for  each  value  of  n  and  a,  correspond¬ 
ing  to  a  discrete  set  of  Y(p=l).  The  results  are  given  in 
table  VT  B  and  in  figures  6.7  to  6.9  below. 
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Third.  Eigensolution  of  (L-2) 
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The  integral  expressions  given  in  table  VIII  B  were  evaluated 
numerically  from 
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APPENDIX  M.  VARIATIONAL  APPROXIMATIONS  TO  ODD  PARITY  NEUTRAL 

PARTICLE  EIGENSOLUTIONS  USING  THE  TRIAL  FJNCTION 

0C  =  R(r)P1(n)+Q(r)P3(n) 


The  comparison  function  -  R(r)P^(n)  +  Q(r)P3(fi) 
reduces  the  Lagrangian  (Vl-la)  to 
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If  R  and  Q,  satisfy  (M-2)  and  (M-3)  then  it  follows,  as  in  the 
case  of  (Appendix  K),  that  satisfies  the  integral 
relations  (V-4)  and  (V-6) .  Note,  however,  that  b^  does  not 
satisfy  (V-7)  since  in  extremising  (M-l)  no  extremisation 
with  respect  to  rectangular  cartesian  scale  parameters  takes 
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place . 


For  large  r,  solutions  of  (M-2)  and  (M-3)  decouple 
and  take  the  form 
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Using  this  information,  solutions  of  (M-2)  and  (M-3)  which 
are  asymptotic  to  zero,  may  be  obtained  numerically  by 
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replacing  the  derivatives  with  central  difference  approximations 
and  solving  the  resulting  difference  equations  as  described 
in  Chapter  V.  The  first  few  eigensolutions ,  obtained  in  this 
way,  are  plotted  in  figures  6.10  to  6.13- 

The  integral  expressions  of  table  VTII  B  take  on 
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First  Eigensolution  of  (M-2)  and  (M-3) 


Third  Eigensolution  of  (M-2)  and  (M— 3 ) 


Fourth  Eigensolution  of  (M-2)  and  (M-3) 
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APPENDIX  N.  VARIATIONAL  APPROXIMATION  TO  ODD  PARITY  EIGENSOLUTION 

WITH  4>d  =  (r2+  br3)  e~arY(p.) 


All  expressions  are  the  same  as  for  (Appendix  L) 
except  that 
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where  B  =  ^  =  cQc2  and  =  cQc^  . 

The  results  for  this  trial  function  are  given  in 
table  VI  B,  while  the  optimal  form  of  Y(d)  is  plotted  in 
figure  6. 14  . 
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APPENDIX  0.  VARIATIONAL  APPROXIMATION  TO  POSSIBLE  ODD  PARITY 
NEUTRAL  PARTICLELIKE  EIGENSOLUTION  USING  THE  TRIAL  FUNCTION 

<t>E  =  R(r)P1(^)  +  Q(r)P3(^i)  +  S(r)P5(|j.) 


Substituting  the  variational  trial  function 
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00 


2  r  1  -n  T  2  ,  1a'2,  1 


M  =  4  /  r  14-R  +  + 


0 


7 


11 


2+  ^(l  +  ^2")  R2+  y(l  +  ^§)Q2 

3  r  r 


+  ^.(i  +  2°)S2-  i  G(R,Q,S)]dr 
r 


(0-1) 


where 


G(R,  Q,S)  =  c-j^rV  c2Q4+  c3S^+  6(c2|R^Q^+  c^S^+c^S^) 


,4 


4 


,2^2 


,2^2.  J2.„2- 


+12 ( c^R^QS  +  CgRQ  S  +  c^RQS  ) 


+  4(c1qR3Q  +  c11RQ"i  +  ci2Q‘3S  +  c13RS-3+  c-^QS0) 


3 


3 


3 


with 


'1 


/  P^U|T  =  ^  , 
0  1 


C2 


,4 


24l 


Jq  P3dn  5Q05  ^ 


1  4 

=  /  P^cljj,  =  .0218450  , 


3 


o 


■  i  j-aht  sh:  :7t:  t  ;  ; ixu;  x  ams.mTFiA'i  jAHTuai: 

UlgiiTla^  .(ti)gH(i)g  +  (ju) ri(/i)a  °,g<t  ■ 


.  • 


*  ,,  ■  -  +  U),  (".}  ■  (  |)XJ.(.?)H  =  -g* 


0» 


0 


1 

} 


‘v 


/ 


3*I9C'v  jS 


> 


. 

•>. 


IV  s  V,,  v 


rW  iw 


0-2 


o 

4^ 

II 

1  p 

/  Pi  p 
0 

2  du  23 

3  ^  315  ' 

°5  = 

1  2 

J  P  P 
0 

^  dq 

0458430  , 

II 

VO 

o 

1  2 

/  Pf  P 
0  0 

c,  dq  =  ■ 

.  0254451  , 

II 

r- 

o 

1  2 

/  P?  F 
0  1 

’  P 

3 

3 

=  .0288600  , 

o 

00 

II 

1 

/  Pn  E 
0 

>2  p 

3 

5  dM- 

=  .0299700  , 

o 

II 

1 

/  Px  E 
0  1 

5  3  F 

'1  dq 

=  .0179820  , 

o 

1 — 1 
o 

1  2 

=  /  pi 

0  1 

P3 

du  = 

2 

35  ' 

1 — 1 

1 — 1 

O 

1 

=  /  Pi 
0 

P3 

3 

dq  = 

12 

3®5  ’ 

OJ 

1 — 1 
o 

1  2 
=  /  p? 
0  3 

P5 

dq  = 

.0199800  , 

QO 

i — 1 

o 

1 

=  J  P1 
0  1 

P1 

dq  = 

.0123406  , 

-=*■ 

1 — 1 
o 

1 

=  /  P3 
0  ^ 

P5 

dq  = 

.0140417  • 

p  (q)  are  the  usual  Legendre  polynomials. 
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The  Euler  -  Lagrange  equations  which  must  be  satisfied 
by  R,  Q  and  S  in  order  that  (O-l)  be  extremised  are: 
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where  GR  =  and  similarly  for  Q  and  S. 

The  resulting  optimised  trial  functions  ^  satisfy 
identically  the  integral  relations  (V-4)  and  (V-6),  for  the 
same  reason  as  the  optimised  trial  functions  and  d>c  (see 
Appendix  K) .  Thus,  the  expression  (O-l)  for  the  mass  may  be 
written  as 

M  =  /  0e  d3x  =  ^  /oG(R^^s)r2dr* 

The  above  set  of  coupled  non-linear  ordinary 
differential  equations  (0-2)  can  be  solved  numerically  using 
the  finite  difference  method  described  in  Chapter  V.  Note  that 
when  r  is  large  the  solutions  of  (0-2)  which  are  asymptotic 
to  zero  are  of  the  form 
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respectively,  while  for  r  «  1  R,  Q  and  S  are  proportional 

•3  R 

to  r,  rJ  and  r^  respectively. 

Numerical  solution  of  the  equations  (0-2)  indicates 

the  existence  of  a  discrete  set  of  orthogonal  (in  the  sense  of 

(VI-l))  eigensolutions  asymptotic  to  zero,  corresponding  to 

a  discrete  set  of  the  constants  A,  B  and  C.  The  resulting  set 

of  optimised  trial  eigenfunctions  then  satisfy  also  the 

orthogonality  relations  (VI-l) :  For  suppose  (R1,Qi,Si)  and 

( R  . , Q . , S . )  are  two  solutions  of  (0-2).  Then  multiplying  the 
J  d  d 

set  of  equations  (0-2)  for  (R^, )  by  Rj,  Qj  and  Sj 
respectively,  performing  the  same  calculation  with  i  and  j 
interchanged,  subtracting  the  resulting  integral  relations 
and  using  Greenes  theorem  we  obtain 

/  Gr  R.  d3x  =  /  gr.ri  d3x 

i  J 

and  similar  results  for  Q  and  S.  These  relations  then  combine 
to  yield 

/  ^  d3x  =  /  d3x  , 

which  is  just  the  relation  (VI-l)  for  the  optimised  trial 
eigenfunctions  bR^ ^ . 
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The  results  for  this  trial  function  are  given  in 
tables  VI -B  and  VIII-B,  while  the  optimal  forms  of  the 
(nodeless)  functions  R,  Q,  and  S  are  plotted  in  figure  6.15* 
The  various  integrals  (p.  95)  which  are  listed  in  table 
VIII -B,  were  evaluated  numerically  from: 
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Nodeless  Elgensolution  of  (0-2) 
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APPENDIX  P.  SOLUTION  OF  FINITE  DIFFERENCE  EQUATION  (V-13) 


In  the  table  below  we  list  a  result  of  iterative 
solution  of  the  finite  difference  approximation  (V-13)  to  (III-4). 
0(r,p)  =  Ak1(r)P1(ix)  for  r  >  5(|^),  with  A  =  10  (|^-)  .  The  mesh 
size  used  is  br  =  .0l(£-)  and  bp  =  .05  •  4>(r,p)  is  plotted,  as 

ill  O 

a  function  of  p,  for  some  (arbitrarily  chosen)  fixed  values 
of  r  in  figures  6.15  and  6.l6  . 
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1  Decreasing  the  mesh  size  (especially-  bp)  had  no  significant 
effect  on  the  results  for  r  >  1,  but  enhanced  instability1  for 
small  values  or  r. 
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APPENDIX  Q.  INTEGRALS  OF  TRIAL  FUNCTIONS  FOR  SPHERICALLY 

SYMMETRIC  CHARGED  PARTICLELIKE  SOLUTIONS 

1.  ^  1  -  xTi 
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H(3)  =  c6(iu3-  |u2+  -  4) 
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G  =  4  c 3  ( -j=u  -  -?j-u  +  -^ ) 
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Gk  =  Z0  O  H(n)  ' 

n=0 


(k)  being  the  usual  binomial  coeficients. 


3.  f(x,Tj,g)  =  1  -  2x°+  r\x2° 
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4.  f(x,T])  =  -i - ii 

x  +1 

c  =  (i  -  1)* 

o  2 
v  =  2cr| 

Defining  H(l)  =  tan 

and  H(n)  =  2(n-l)  ^  CTin~1+  ( 2n  "  3)  H(n-l)]  when  n>l, 

F(0)  =  j  c3 
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G,  =  2  (~l)n  (k)  F(k-n)iin,  where  (j*)  are  binomial 
k  n=0  n 

coeficients . 

5.  f(x,Ti,g)  =  cos(x+ri)  +  a 

i  7T  .  -1_ 

s  =  c  +  t)=2‘  +  sin  o 

V  =  (1  -  o-^)2 

1  -|  o  -j  ■  13121 
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APPENDIX  R.  APPROXIMATE  ‘SOLUTION*  OF  ( VII -13) 


Figures  7.7  and  7.8  below  are  plots  of  cross-sections 
(for  certain*  arbitrarily  chosen*  fixed  values  of  r)  of  a 
numerically  obtained  "solution"of  (VII-13) •  The  method  used  in 
obtaining  these  results  is  the  method  of  iteratively  solving 
the  finite  difference  approximation  to  ( VII -13);  this  method 
is  described  in  Chapter  V. 
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me 

For  the  example  plotted  below  a  =  -2.177  — — , 

2  & 
b  =  2.77  —  and  p  =  .1043  %:•  The  mesh  size  used  is 

g  o111 

5r  =  ,oi  ~  and  6|jl  =  .05.  The  results  were  found  to  be 
essentially  unaltered  when  the  mesh  size  was  varied  (decreased) 
except  for  the  region  where  r  <  .4*  where  the  unstable 
"oscillations "which  are  apparent  (fig.  7-8)  grew  as  the  mesh 
size  was  decreased  (especially  6|a)  . 
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